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BRIEF  OUTLINE  OF  RESEARCH  FINDINGS 


This  summary  describes  the  work  on  ARO  Grant  Number  DAAG29-76-G- 
0331  from  24  September  1976  to  23  September  1979.  The  purpose  of  this 
grant  is  to  develop  the  basic  theory,  algorithms  and  computational  tech¬ 
niques  to  analyze  the  electrically  thin,  dielectric  loaded,  cavity  backed 
radiator,  commonly  called  a  microstrip  antenna.  The  overall  approach 
will  be  to  develop  techniques  which  are  applicable  to  a  wide  class  of 
geometries,  and  to  apply  these  general  techniques  to  this  specific  antenna 

Our  approach  can  be  divided  into  six  steps  as  follows: 

1.  Develop  a  general  purpose  moment  method  code  for  analyzing 
geometries  involving  wires,  plates  and  wire  plate  attachments. 

2.  Show  that  this  general  purpose  code  is  applicable  to  the 
electrically  thin  cavity  antenna,  or  microstrip  antenna. 

3.  Adapt  the  technique  to  treat  patches  of  "arbitrary"  shape. 

4.  Adapt  the  technique  to  treat  attachments  near  edges. 

5.  Adapt  the  technique  to  treat  dielectric  loading. 

6.  Verify  accuracy  by  comparison  with  measurements. 

At  present  all  steps  except  number  4  are  complete.  Thus,  we  are 
able  to  analyze  coax  fed  microstrip  antennas  of  rectangular  and  non- 
rectangular  shape.  Although  considerable  effort  was  applied  to  the  study 
of  wires  attached  to  edges  of  plates,  we  do  not  have  a  satisfactory 
technique  for  treating  edge  fed  microstrips.  Our  current  techniques 
predict  resonant  frequency,  but  not  impedance  level  accurately  for  the 
edge  fed  microstrips. 
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The  ipp^niices  reproduce  our  journal  publications  resulting  from 
this  research.  The  application  of  these  techniques  to  the  microstrip 
antenna  is  the  subject  of  a  manuscript  now  in  preparation.  It  is  sum¬ 
marized  below. 

Most  previous  work  on  microstrip  antennas  can  be  divided  into  three 
categories.  First,  there  are  the  transmission  line  models,  where  the 
unknowns  are  the  propagation  constant,  characteri Stic  impedance,  and 
loads  on  an  equivalent  transmission  line  which  approximates  the  microstrip 
antenna.  Second,  there  are  the  modal  solutions,  where  the  unknowns  are 
the  distribution  of  cavity  modes  in  the  cavity  formed  by  the  region  between 
the  microstrip  patch  and  the  ground  plane.  By  contrast,  in  our  solution 
the  unknown  is  the  current  flowing  on  the  microstrip  patch.  This  current 
is  determined  by  using  the  method  of  moments  to  solve  an  integral  equation 
in  which  the  current  is  the  unknown. 

The  basic  building  block  in  our  solution  is  a  moment-method  surface- 
patch  modeling  technique  for  analyzing  bodies  composed  of  wires  and  con¬ 
ducting  plates.  First  it  will  be  shown  how  the  basic  surface-patch 
technique  is  applied  to  the  air-dielectric  microstrip  antenna  of  Figure 
la.  The  antenna  consists  of  a  plate  or  patch  a  distance  T  above  a  ground 
plane.  The  antenna  is  coaxially  fed  with  the  center  conductor  extending 
beyond  the  ground  plane  and  contacting  the  top  patch.  In  Figure  lb  the 
ground  plane  is  removed  and  the  images  of  the  top  patch  and  the  center 
conductor  are  inserted.  The  coax  aperture  is  modeled  by  a  voltage  source. 
Previous  experience  by  this  author  and  others  indicates  that  the  coax 
center  conductor  need  not  be  modeled  as  a  metallic  wire,  but  can  be  modeled 
rather  simply  as  a  constant  current  filament  where  for  convenience  a 
current  of  one  amp  is  used.  This  is  shown  in  Figure  lc.  The  geoemtry 
of  Figure  lc  can  be  analyzed  with  a  surface-patch  moment-method  code. 
Briefly,  the  current  on  the  patches  is  expanded  in  terms  of  N  expansion 
modes: 
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Figure  1.  Models  for  the  microstrip  antenna. 


Typical  rectangular  surface  patch  current  expansion  or  testing 
modes  are  shown  in  Figure  2a.  The  moment-method  solution  results  in 
a  system  of  simultaneous  linear  equations  of  the  form 

f zl  I  =  V  (2) 

4.  J 

r  i 

where  IZj  is  the  NxN  impedance  matrix,  I  is  a  column  vector  containing 
the  N  expansion  coefficients  in  Equation  (1)  and  V  is  the  excitation 
column  vector.  A  typical  element  in  is  given  by 

zmn  *  -  1  in  '  f  (3) 

n 

where  Em  is  the  free  space  fields  of  the  mth  testing  function  and  the 

^  h 

integration  is  over  the  surface  of  the  nLr  expansion  mode.  A  typical 
element  of  V  is  given  by 


Vm  ’  f  ii  *  ^  ds 

i 


(4) 


where  is  the  impressed  one-amp  current  in  Figure  lc.  Application 
of  this  technique  is  in  principle  straight-f orward,  however,  in  practice 
it  is  challenging  since  certain  terms  in  the  moment  method  impedance 
matrix  need  to  be  evaluated  to  within  about  0.1%  accuracy. 


Mext  it  will  be  shown  how  the  moment  method  solution  is  modified 
by  the  dielectric  slab  between  the  microstrip  patch  and  the  ground  plane. 
Figure  3a  shows  the  model  of  Figure  lc,  but  with  the  addition  of  the  di¬ 
electric  slab  and  its  image.  Using  the  volume  equivalence  theorem,  the 
slab  may  be  removed  and  replaced  by  the  equivalent  volume  polarization 
currents 


iv  =  >(e-:0)E  (5) 

where  :  is  the  permittivity  of  the  slab,  cQ  is  the  permittivity  of  free 
space,  and  £  is  tne  actual  electric  field  intensity  in  the  slab.  This 
is  illustrated  in  Figure  3b.  Jy  is,  of  course,  unknown  since  E^  is  unknown. 
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Figure  2.  Surface  patch  expansion  and  test  modes. 


Microstrip  antenna  with  dielectric  slab 


We  have  had  previous  experience  in  modifying  mjment-method  solutions 
to  account  for  dielectric  inhomogeneities  where  iJ  and  Js  are  considered 
either  as  independent  unknowns  or  as  dependent  unknowns. 


Here  we  consider  them  to  be  dependent  unknowns.  In  this  case 
Equation  (2)  is  modified  to 


fz+Azll  =  V  (6) 


where  AZ  is  a  matrix  which  accounts  for  the  presence  of  the  dielectric 
slab.  A  typical  term  in  [*zl  is  given  by 


=  -  f  in  •  r  dv 


where  E_^  is  the  field  of  the  n*'1'1  expansion  mode  in  the  presence  of  the 
dielectric  slab,  Em  is  still  the  free  space  field  of  the  m^11  test  mode, 
and  the  integral  is  over  the  volume  of  the  dielectric  slab. 


If  the  exact  Green's  function  for  the  slab  is  used  in  evaluating 
E_",  then  the  presence  of  the  slab  is  handled  exactly.  Unfortunately, 
the  Green's  function  for  the  slab  involves  the  Sommerfeld  integrals  which 
are  difficult  and  time  consuming  to  evaluate.  Thus,  to  have  a  numerically 
efficient  technique  for  evaluating  the  ^Z  matrix  we  need  a  much  simpler 
method  for  evaluating  E*].  The  first  approximation  used  was  that  only  the 
z  the  component  of  E^J  would  be  significant  in  the  slab.  Next,  previous 
work  suggested  that  z  •  E^  in  the  slab  could  be  reasonably  approximated 
by  the  z  component  of  the  fields  radiated  by  in  a  homogeneous  medium 
of  permittivity  c.  This  was  found  to  be  a  reasonable  approximation 
except  that  it  ignores  the  fact  that,  for  a  surface  current  on  the  surface 
of  a  slab,  the  normal  component  of  the  D  field  in  air  will  be  somewhat 
less  than  the  normal  component  of  the  D  field  in  the  slab.  This  effect 
was  taken  into  account  using  the  method  of  dielectric  images  and  considering 
the  fields  to  be  quasi-static  in  the  slab. 
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Figure  A  shows  a  comparison  of  measured  and  calculated  input  im¬ 
pedance  for  a  coaxial ly-fed  rectangular  microstrip.  Although  the 
details  were  not  mentioned  above,  losses  in  the  dielectric  and  in  the 
finitely  conducting  patches  were  included  in  the  calculations,  and  have 
a  significant  effect  on  the  bandwidth  and  the  impedance  level.  Non- 
rectangular  plates  can  be  treated  using  the  modes  of  Figure  2b. 
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Figure  4.  Input  impedance  of  rectangular  microstrip  antenna. 


APPENDIX 


Rigorous  Near  /one  Field  E  xprexsions  fur  Rec  tangnlar 
Sinusoidal  Surface  Monopole 

iV  K  II  KI<  HVIOND.  SENIOR  VDMHER.  IE.EE.  U  VI  POZAR 
AND  E  OW \RD  H  NO  VI  AS  VI E  VI  HE  R.  IHI 


\h\tran  Ihe  near  rone  fields  of  che  sinusoidal  surface  monopolc 
are  prevnted  in  lerms  of  exponential  integrals. 


I  INTRODUCTION 

It  is  widely  known  that  a  sinusoidal  electric  line  source  with 
arbitrary  length  has  simple  rigorous  near-zone  field  expressions 
(  I  |  f  urthermore  this  has  been  used  with  considerable  advan¬ 
tage  in  moment-method  solutions  for  thin-wire  antennas  and 
scatlerers  with  arbitrary  shape  [2|-[4|  In  moment-method 
solutions  for  sheet-metal  antennas  and  scatterers.  it  is  advan¬ 
tageous  to  use  surface  dipoles  (instead  of  thin-wire  dipoles)  for 
expanding  the  unknown  surface  current  density  in  a  basis  set 
|  5  |  A  great  deal  of  computational  expense  can  be  avoided  by 
selecting  a  surface  current  mode  with  near-zone  fields  expres¬ 
sible  as  a  finite  series  of  tabulated  functions.  Only  the  sinus¬ 
oidal  surface  current  distribution  is  known  to  have  this 
important  property  The  field  of  the  sinusoidal  dipole  may  be 
obtained  by  summing  the  iieius  of  iwo  sinusoidal  monopoies. 

Using  the  field  expressions  peresented  here,  computer  sub¬ 
routines  have  been  developed  for  the  fields  of  sinusoidal  sur¬ 
face  monopoles  and  dipoles.  These  algorithms  are  found  to  be 
efficient  and  stable  even  when  the  observation  point  is  near  or 
on  the  source  Furthermore,  the  calculated  fields  satisfy  all  the 
appropriate  boundary  conditions  at  the  source.  When  the 
distance  to  the  observation  point  is  not  too  small,  the  fields 
may  be  obtained  readily  with  numerical  integration.  In 
moment-method  applications,  however,  the  observer  is  often 
on  or  extremely  close  to  the  source  where  the  numerical 
approach  is  difficult. 

II.  NEAR  ZONE  FIELD  EXPRESSIONS 

C  onsider  a  time-harmonic  source  with  the  following  electric 
surface  current  density: 


( / »  sinh  y(z2  —  z)  +  /2  sinh  y(z  —  Zj  )|  cosh  yy 

■f  *  - 

2  sinh  iyh)  sinh  (yw) 


Ihe  time  dependence  i,'~'  is  understood  As  illustrated  in  Fig. 
1.  this  source  lies  in  the  i:  plane  in  the  Tegion  »  ■-  t  <"  w 

and  .'|  .'  ^  I  is  a  planar  rectangular  source  with  edge 

dimensions  and  h.  where  h  ~  ;2  Cj ,  The  constants 
and  /2  denote  the  terminal  currents  in  amperes  at  and  z2. 
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I  ig  1  Geometry  of  rectangular  surface  inonopule  and  observation 

point 

respectively.  Let  this  rectangular  surface  monopole  radiate  in  a 
homogeneous  conducting  medium  with  intrinsic  impedance  rj 
and  propagation  constant  y  By  integrating  (with  respect  to 
v  )  the  fields  of  the  sinusoidal  line  source,  we  obtain  the 
following  expressions  for  the  fields  of  the  sinusoidal  surface 
monopole: 

tx(x.  y.  ;)  =  jr\C[{l2S2  -  /,Sj )  sinh  yh  +  /j  2E,  *  /21£2  I 

C> 

Ey(x,  y,  z)  =  r)C((/2  T2  ~  I  \TX)  sinh  y/i  +  /j  20'i  +  /2  ,o'2  j 

1 3 ) 


£,(x.  y.  ; 

: )  —  (Ei  2  *  l  X2  j  E  2 ) 

i  4  i 

Hx(x.  i. 

z )  =  C[(/2G2  ~  / 1 G i  )  sinh  y h  +  li2£i 

*  f2 1  E 2  ) 

<5> 

Hy(x.  y. 

z )  =  /Clf/jEj  —  /2E2  1  sinh  yh  -/j2.S'i 

/2  1-^2  1 

1 6 » 

Hz{  x.  y. 

z  1  =  0 

( 7 1 

S~~ 

II 

2ti  sinh  yh  sinh  yu  ) 

i  8) 

fi  2  ~  h 

cosh  yh  —  /2 

(91 

/2 j  =  l2 

cosh  yh  —  / j 

(10) 

=  2 

klmt(k.  1,  m,  ti) 

till 

"n  2 

>»!  E  ( k.  1.  >11.  II ) 

1 1:> 

2 

k.'E [k.  1.  in.  n | 

i  1 3  i 

£„  -  3E(0.  1 . 0.  it)  -El 0.  1 .  0.  n)  3  2^-*^  f  *0 


(14) 

y„  =  e.(o  i . 

0.  n)  -k(0.  -1.  0.  n) 

(15) 

£(k.  1.  m.  <t) 

=  e,*>V,',ve'n><I  /  _ 

1  lb) 

J. 
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I  ig  2  Geometry  of  rectangular  surface  dipole  with  fields  given  in 
Table  I 


TABLE  I 

.  t  PCTRIC  FIELD  STHLNGTH  IN  V  M  fOR  SURf ACE  DIPOLE 
IN  HO.  2  ALONG  LINE  (\,  y,  i )  =  lx.OOSX.  0  1M 


*  =  Rn  1  +  /A.  A  7-  l(y  -  y\  )  +  m(z  ~  Zn  ) 
6  =  R n2  +  lkx  +  l(y  ~  >2)  +  <*(Z  ~  Zn) 


1 

2 


i-~  1 


(17) 

(18) 

(19) 


As  indicated  in  Fig.  1 .  R„ ,  and  Rn2  denote  the  distance  from 
a  corner  of  the  source  to  the  observation  point  In  (17)  and 
(18).  i-,  =  u  and  i2  -  ■*  For  the  exponential  integral  in 
I  16),  subroutine  F.XPJ  isavailable  [2],  (3],  In  the  summations, 
the  indices  Jc.  /,  and  m  take  on  only  the  values  -1  and  +l,so 
there  are  eight  terms  altogether  in  the  triple  sum.  The  expres¬ 
sion  for  Ey  -is  correct  only  if  x  differs  from  zero.  When  x  van¬ 
ishes,  however,  the  correct  field  may  still  be  obtained  from  (3) 
with  any  suitably  small  value  of  x  such  as  x  =  w/106 

The  sinusoidal  monopole  has  line  charges  on  the  edges  at 
:  C(  and  c  ~  .'2.  but  they  will  disappear  when  another 

monopole  is  connected  to  one  edge  I  to  form  a  sinusoidal 
dipole)  and  the  current  al  [he  other  edge  is  set  equal  to  zero 
For  this  reason,  the  contributions  from  the  line  charges  are 
omitted  in  the  field  expressions. 

Fig.  2  illustrates  a  surface  dipole  in  the  y:  plane  with 
terminals  at  z  =  0  and  one  ampere  terminal  current.  At  300 
MHz,  the  electric  field  intensity  of  this  dipole  (in  free  space)  is 
listed  in  Table  I. 
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APPENDIX  3 


Electromagnetic  Modeling  of  Composite  Wire  and  Surface  Geometries 

[  H  NFWMAN,  MFMBFK,  IFFF  \NI)  I)  M  »'()/  \  R.  si  I  l)F  M  MhMKIK.  Ithf 


4t>aract  A  moment  method  solution  to  the  problem  of  radiation 
or  scattering  from  geometries  consisting  of  open  »t  closed  surfaces, 
wires,  and  wire  surface  junctions  is  presented,  rhe  method  is  based  on 
the  sinusoidal  reaction  formulation.  Several  examples  of  input  imped 
ance  calculations  illustrate  the  versatility  .  accuracy,  and  computational 
efficiency  of  the  method. 

I  INTRODUCTION 

HI-  MF.THOD  of  moments  is  perhaps  the  most  widely  used 
tool  for  the  electromagnetic  modeling  of  bodies  which  are 
not  large  in  terms  of  a  wavelength  The  advantages  of  moment 
methods  are  accuracy,  versatility,  and  the  ability  to  compute 
near-  as  well  as  far  cone  parameters. 

The  most  widely  used  forms  of  the  method  of  moments  are 
the  thin-wire  computer  programs  [  1  ]  — ( 3 1 .  These  programs  are 
ideal  for  modeling  most  wire  antennas  By  forming  a  wire-grid 
model  they  also  can  be  used  to  model  solid  surfaces  Unfortu¬ 
nately,  the  size  of  the  surface  which  can  be  modeled  practically 
is  severely  limited  since  it  requires  many  wires  to  accurately 
model  a  solid  surface.  Further,  even  a  very  fine  wire  mesh  may 
not  yield  accurate  near-zone  parameters,  such  as  current  dis¬ 
tribution  or  impedance 

More  recently,  surface-patch  models  have  been  developed 
for  modeling  solid  surfaces  The  advantage  of  the  surface  patch 
solution  is  that  fewer  unknowns  are  required  per  square  wave¬ 
length  of  surface  area  Wang  et  al  |4)  have  developed  a  model 
for  wires  and  plates  based  on  the  s.nusoidal  reaction  formula¬ 
tion  Alhertven  i  r  j/  (5)  used  pulse  expansion  modes  to  model 
wires  and  plates,  including  the  case  where  a  wire  touches  a 
plate  Burke  and  Poggio  have  incorporated  these  results  into  a 
user-oriented  computer  code  [3]  Parhami  et  ill.  have  treated 
the  problem  of  a  wire  surface  junction  using  the  finite  differ¬ 
ence  technique  [61 

The  purpose  of  this  paper  is  to  present  a  technique  for 
modeling  wires  and  surfaces,  including  the  case  where  the 
wires  contact  the  surfaces,  which  is  sufficiently  accurate  to 
compute  impedance.  The  technique  is  based  on  the  sinusoidal 
reaction  formulation  and  is  applicable  to  open  as  well  as  closed 
surfaces  Continuity  of  current  is  enforced  on  the  wires,  on  the 
surfaces  and  at  the  wiresurface  junctions.  The  accuracy  and 
versatility  of  the  model  is  demonstrated  by  several  numerical 
examples  Attachment  points  are  restricted  to  be  somewhat 
removed  from  an  edge 

II  THEORY 

Pie  Reaction  WethoJ 

Consider  a  geometry  of  arbitrarily  shaped  scatterers  in  a 
homogeneous  medium  as  shown  in  Fig.  1 .  S  is  the  surface 
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ScatTERgus 

Fig  1  General  problem  geometry 


enclosing  the  scatterers.  and  n  is  the  unit  outward  normal  to 
S  The  scatterers  may  consist  of  solid  surfaces  and  thin  w;res 
and  the  wires  may  contact  the  surfaces  For  now  the  surfaces 
will  be  considered  closed.  The  sources  (I,,  M,  i  generate  the 
fields  (E,  H)  in  the  presence  of  the  scatterers.  The  eJUJ'  time 
dependence  is  suppressed  From  the  surface-equivalence 
theorem  the  field  interior  to  the  surface  S  will  vanish  without 
changing  the  extenor  fields  if  the  surface  current  densities 


J,  =  h  X  H, 

U) 

M,  =  E  x  A 

1 

are  introduced  on  the  surface  5  The  scattering  obstacles  may 
then  be  replaced  by  the  ambient  medium  without  altering  the 
field  anywhere.  The  scattered  field,  radiated  by  t  Js .  Ms  1  in  the 
ambient  medium,  is  defined  as 

E.  =  E  —  E, 

(3) 

H,  =  H  — H, 

(4) 

An  electric  test  source,  which  radiates  the  fields  iEj-.  Hri 
in  the  ambient  medium,  is  now  placed  in  the  interior  region  of 
S  Noting  that  this  test  source  has  zero  reaction  with  the  sources 
(J,.  M,  1  and  iJ,.  M, ).  the  reaction  integral  equation  [4]  is 
obtained 

f  jO,  ■  Er  -  M.  •  Ht)  ds  +  If  fo,  •  Et  —  M,  •  Hr)  Jv  =  0 

S  V 

1 5 ) 

where  the  volume  integral  is  over  the  source  volume  The  reac¬ 
tion  equation  1 5  Ms  used  to  determine  the  unknown  surface 
currents  (  K.  M,  l  For  conductors  of  finite  conductivity  the 
impedance  boundary  condition 

M,  =  Z,i,  X  A  (6) 

is  used  where  Z ,  may  be  a  function  of  position  For  simplicity, 
only  perfect  conductors  are  considered  here,  thus  /,  =  M,  =  0 
In  deriving  (5)  it  was  assumed  that  all  surfaces  were  closed 
so  that  Schelkunoffs  equivalence  theorem  could  be  used 
However,  it  can  be  shown  that  (5)  applies  equally  well  for 
open  surfaces.  An  open  surface,  such  as  a  fictitious  plate  of 
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zero  thickness.  an  he  considered  to  he  a  limiting  case  of  a  real 
plate  of  finite  thicknes  as  the  plate  thickness  goes  to  zero  In 
general,  different  currents  exist  on  the  top  and  bottom  surfaces 
of  the  real  plate  As  the  plate  thickness  goes  to  zero,  the  fields 
radiated  by  the  top  and  bottom  currents  become  identical  to 
the  fields  radiated  by  a  single  surface  current  located  on  the 
plate  center  This  single  surface  current  is  the  vector  sum  of 
the  top  and  bottom  surface  currents  and  is  the  current  which 
must  be  determined  to  treat  open  surfaces  If  electric  test 
source.,  are  used  in  ( 5  t.  then  J,  will  be  the  vector  sum  of  the 
current  on  the  top  and  bottom  surfaces  ( 7 J 

The  integral  equation  (5)  is  solved  by  the  moment  method. 
The  unknown  current  J,  is  expanded  in  a  set  of  V  basis  (ex¬ 
pansion  I  functions, 

lnK.  (  7 ) 

n  -  l 

and  i  5  (  is  enforced  for  V  electric  test  sources  placed  in  5.  Thus 
(5 )  reduces  to  the  set  of  simultaneous  linear  equations 

£  tnZmn  =  Vm.  m  =  1,2  "IV  (8) 

n  =  I 


where 


-<i 

- 1, 

-1 - 51 - * - *■  1 

»  V 

(b) 


z«0 

(c) 

Fig  2.  Expansion  and  test  modes,  (a)  Thin-wire  f-dipole.  (b)  Sinus¬ 
oidal  surface  F-dipole  (c)  Attachment  dipole 
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m  n 


ds 


liing  the  thin-wire  approximation  the  current  on  this  dipole 
(9)  is  given  by 


n 

r  r  r 

J."  =  — 

sin  k(z  —  z j )  ^  sin  k(z3  —  z)~ 

Pi  +  P2 

j  J  J  ( J(  •  Em  —  M,  •  Hm  )  dv 

U0) 

2tt  a 

sin  *(z2-z1)  sin  k(z3  ~  z2) 

v 


where  *  Em  HJ  are  the  fields  of  the  "r th  test  source  radiating 
in  the  medium  (p.  e)  and  the  integration  in  (9)  is  over  the 
surface  of  the  nth  expansion  mode 

The  expansion  and  test  modes  used  here  are  identical  Thus 
the  method  is  a  Galerkin  method,  and  a  symmetric  impedance 
matrix  |/|  results  The  form  of  the  expansion/test  functions 
will  now  be  defined  It  is  this  form  which  ultimately  determines 
the  accuracy  and  efficiency  of  the  solution  as  well  as  the  types 
of  geometries  which  may  be  modeled. 


where  Pt  and  P2  represent  pulse  functions  with  unit  value 
when  zt  <  z  <  z2  and  z2  <  z  <  Z3,  respectively,  and  are  zero 
elsewhere.  Also,  a  is  the  wire  radius  and  k  =  2n/\.  This  choice 
of  mode  has  the  advantages  that  the  near-zone  fields  and  the 
7.mn  are  known  in  closed  form  (8]  These  modes  are  placed  in 
an  overlapping  array  on  the  wire  ensunng  continuity  of  current 
on  the  wire 

2}  Surface-Patch  Mode  The  surface-patch  mode  is  a  surface 
I'-dipole  consisting  of  two  sinusoidal  surface  monopoles  A 
surface  F -dipole  with  interior  angle  of  ISO3  is  shown  in  Fig. 
2(b).  The  current  on  this  dipole  is  given  by 


E  xpansion  I Testing  Functions 

Three  basic  types  of  modes  are  used  wire  dipole  modes, 
surface  dipole  modes,  and  a  special  attachment  mode  when¬ 
ever  a  wire  connects  to  a  surface  With  this  choice  of  functions 
geometries  consisting  of  flat  surfaces  thin  w  ires,  jnd  wire-surface 
connections  may  be  modeled  \  piecewise  flat  approximation 
can  be  made  to  model  singly  curved  surfaces  Note  that  all  of 
the  modes  will  involve  sinusoids  with  the  free  space  wave- 
number  This  allows  all  but  one  of  the  integrations  required  to 
find  the  Em  to  be  done  in  closed  form,  thus  avoiding  the  very 
difficult  l/r3  singularity  ot  second  derivative  associated  with 
findirg  the  extreme  near-zone  fields  of  an  electric  current 
source,  as  is  required  to  evaluate  self  or  overlapping  impedance 
elements  The  one  integration  not  available  in  closed 

form  is  associated  with  the  disk  component  of  the  attachment 
mode,  as  discussed  below 

I )  Thin  Wire  Mode  The  wire  mode  used  is  identical  to  that 
used  by  Richmond  [  1 1  It  is  a  piedewise-sinusoidal  F-dipole 
consisting  of  two  sinusoidal  monopoles  A  1-dipole  with  a 
180'  internal  angle  lying  on  the  r-axis  is  shown  in  Fig.  2(al. 


kPx  sin  *(r  —  zt)  cos  ky 

J/  =  z -  +  i 

2  sin  k(z2  ~  Zj)  sin  kw 


kP2  sin  k(z3  —  z)  cos  ky 
2  sin  k(z3  —  z2)  sin  kw 
(12) 


where  l\  and  P,  represent  unit  pulse  (unctions  as  before  This 
mode  is  similar  to  that  used  by  Wang  •  t  al  [4  |  except  that  the 
sinusoidal  variation  is  with  the  speed  of  light  and  a  cosine  vari¬ 
ation  is  used  transverse  to  J ,*  instead  of  the  constant  variation 
used  by  Wang  el  jl  Because  of  this  the  fields  due  to  the  current 
of  (12)  may  be  found  in  closed  form  |9] .  Two  orthogonal  and 
overlapping  arrays  of  the  surface-patch  modes  are  placed  on 
the  surface,  allowing  a  two-dimensional  vector  surface  current 
density.  Results  presented  in  the  next  section  indicate  that 
reasonable  accuracy  can  be  obtained  with  2w  or  :2  :2  or 

z3  z2  as  large  as  0  25  X. 

3)  Attachment  Mode  When  a  wire  is  attached  to  a  surface 
a  special  attachment  mode  is  introduced.  The  purpose  of  the 
attachment  mode  is  twofold  first  is  to  establish  continuity  of 
current  at  the  wire/surface  junction  and  second  is  to  insure 
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rha:  in  the  immediate  vicinity  of  the  attachment  the  surface 
current  density  has  the  proper  p  polarization  and  1  p  depend 
ence  Hus  mode  consists  of  two  parts  a  wire  monopole  and  a 
disk  monopole,  as  shown  in  Fig  2(c).  The  wire  monopole 
cuirent  density  is  similar  to  the  thin-wire  monopole  mode 

1  sin  k(z,  —  z) 

j  w  _ -  -  *  i 

•  7  ’  (13) 

-tra  sin  *z2 

Thus 

sin  -  p) 

•*.  P,  a<p<.b  (14) 

2trp  sin  k(b  —  a) 

where  a,  b  are  'he  mner  and  outer  radii  of  the  annulus  Note 
that  the  total  current  on  the  disk  at  p  =  a  is  equal  to  the  total 
current  on  the  wire  segment  at  z  =  0  insuring  continuity  of 
current  at  the  attachment  Also  observe  that  the  sin  kib  p ) 
function  in  the  numerator  of  ( 14)  forces  the  disk  current  to  be 
zero  at  the  disk  edge  p  =  h  It  is  this  property  which  allows 

the  disk  to  be  placed  on  the  surface  and  still  maintain  conti¬ 

nuity  of  current  on  the  surface  The  sin  kib  p)  function  was 
chosen  (rather  than,  say,  *  p  or  smart  6  p ))  since  this  per¬ 
mits  the  fields  of  the  disk  to  be  obtained  with  only  one  numer¬ 
ical  integration 

The  attachment  mode  is  applied  by  placing  it  directly  on 
the  surface-patch  modes  uherever  the  wire  meets  the  surface 
Attachment  points  are  not  restricted  to  he  in  the  center  of  or 
at  the  corners  of  the  ,urtace  patch  modes  However,  a  more 
detailed  treatment  of  the  attachment  mode,  including  the  edge 
singularity,  would  be  required  tor  attachment  points  less  than 
about  0  1  A  from  an  edge  The  inner  radius  a  of  the  annulus 
corresponds  to  the  wire  radius,  and  results  presented  below 
indicate  that  the  outer  radius  b  has  little  effect  on  the  final 
result  if  it  is  chosen  to  be  between  0.1  A  to  0. 25  A  (see  Fig.  3 ) 

With  these  expansion  functions,  current  continuity  in  the 
direction  of  J,  is  always  maintained,  although  the  current  may 
not  he  continuous  in  directions  orthogonal  to  the  direction  of 
current  flow  Thus  7  •  J,  =  --/uiPj  is  always  finite,  and  no  line 
charges  appear  on  the  surface 

For  antennas,  feeds  may  be  placed  at  the  endpoints  of  any 
wire  segment,  including  the  attachment  segment.  For  example, 
a  monopole-on  a  ground  plane  may  be  modeled  by  inserting  a 
generator  at  the  base  of  the  monopole  between  the  attachment 
mode  wire  segment  and  disk  In  the  work  presented  here  the 
delta-gap  feed  was  used  although  the  magnetic  frill  model  also 
could  he  used 

The  scatterer  or  antenna  geometry  is  thus  modeled  by  Vu 
wire  modes,  V,  surface-patch  modes,  and  V4  attachment 
modes  Then,  V  Vu  +■  V,  +  V4,  and  <8 1  is  solved  for 
the  unknown  current  samples  /„  by  standard  Gauss-Jordan 
elimination 

III  NUMERICAL  EXAMPLES 

Several  numerical  examples  will  now  be  presented  which 
illustrate  the  accuracy  and  versatility  of  the  solution  presented 
in  Section  II  In  each  case,  the  number  of  modes  used  is  sum¬ 
marized  in  Table  1  Note  that  N,  is  the  total  number  of  surface- 
patch  dipoles.  Thus  there  are  ;V,/2  dipoles  in  each  of  the  two 
orthogonal  directions  on  the  surface  The  data  presented  are 
input  impedance  or  admittance,  since  this  tends  to  be  a  sensi¬ 
tive  indicator  of  overall  accuracy 


o  oi  a  z  os 

disk  saoius  ^4 


Fig  )  Input  admittance  lor  geometry  of  Fig  4  with  \  =14  versus 
disk  radius  /  -  160  MHz 

TABLE  1 

P  SRAM)  TfcKS  OF  CALl_l  LAIIONS 


f  ISUBE  IJIVJ  X 

to 

j  %  '  <|ttta  .wtxovus 

3/  -  T' 

Kz  f".  -  -s 

4/ -  2 

40  l  43 

4/  -  1 

12  4 

*/ - 1  IZI 

121 

5/ .  2 

34  J  3? 

6/  •  •  •  •  •  2 

21  24 

9/  -  4 

12  2  18 

9/  -  7 

77  !  3  3? 
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In  the  experimental  models  the  wires  were  brass  rods. 
0.0008  m  in  radius,  with  a  conductivity  of  15  X  10®  mho  m 
The  finite  conductivity  of  the  wires  was  taken  into  account  in 
the  manner  used  by  Richmond  t  \ 01 .  The  metal  plates  were 
made  of  brass  or  aluminum,  and  their  finite  conductivity  was 
not  taken  into  account  In  all  cases  the  driven  wire  element 
was  driven  with  a  coaxial  feed  with  outer  radius  0  004’  m  All 
measurements  were  made  by  the  authors 

The  calculated  results  were  all  made  using  the  same  general 
purpose  computer  code  No  advantage  of  inherent  symmetries  in 
the  geometries  was  made.  The  resulting  symmetries  in  the  expan¬ 
sion  mode  currents  offered  a  useful  check  on  the  calculations 
Eig.  4  presents  the  input  admittance  for  a  monopole  on  a 
square  finite  ground  plane  versus  frequency  in  the  vicinity  of 
the  first  resonance  The  calculation  was  done  with  V  =  14  and 
with  43  unknowns  and  is  compared  with  measured  data  and 
with  a  wire-gnd  model  which  used  121  unknowns  This  example 
shows  the  accuracy  and  rapid  convergence  of  the  pat-h  io.a- 
non  as  well  as  a  savings  in  storage  over  the  wire  grid  "■.ode. 
The  run  time  for  the  N  =14  case  »ii  about  2;  min  n  a 
Datacraft  6024  3,  which  is  equivalent  to  about  2-3  min  on  an 
IBM  3'0-165.’  In  comparison  the  wire-grid  model  took  about 
40  min  on  the  Datacraft  6024  3 

Fig  3  shows  the  effect  of  the  attachment  disk  radius  b  on 
the  input  admittance  of  the  same  geometry  as  Fv*  4  This 
graph  shows  that  the  disk  radius  is  not  cntical  and  usually  can 
be  selected  between  0  I  A  -»  0.25  A 

The  impedance  versus  height  of  a  monopole  on  a  1  4  X 
square  ground  plane  is  shown  in  Fig  5  These  calculations  are 
compared  with  the  input  impedance  of  a  monopole  of  the 
same  height  on  an  infinite  ground  plane  computed  by  image 

1  Recent  improvements  in  the  technique  and  computer  code  have 
reduced  the  run  time  by  about  a  factor  of  6.  while  increasing  accuracy 
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Fig-  4  Input  admittance  versus  frequency  for  monopole  (ft  =  0.421  m. 
a  =  0  0008  m)  on  square  ground  plane  <9.914  x  0.914  ml,  com¬ 
pared  with  measured  data  and  wire-gnd  modeling  b  -  0.2  X. 


Fig  5  (a)  Input  resistance  versus  monopole  height  for  monopole 

centered  on  1  4  X  square  ground  plane  compared  with  same 
monopole  on  infinite  ground  plane  as  computed  by  image  theory 
lb)  Input  reactance  versus  monopole  height  for  monopole  centered 
on  1  4  X  square  gound  plane  compared  with  same  monopole  on 
infinite  ground  plane,  as  computed  by  image  theory 
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hig  t>  Input  admittance  versus  frequency  for  geometry  of  Fig  4 
with  attached  reflected  plate  compared  with  measured  values 

theory  Since  the  I  4  X  ground  plane  is  fairly  large,  good  agree¬ 
ment  is  obtained 

Fig  6  shows  the  input  admittance  for  a  geometry  identical 
to  that  of  Fig.  3,  except  that  a  reflecting  plate  is  attached  to 
one  side  of  the  ground  plane.  Note  that  the  width  of  the  re¬ 
flecting  plate  differs  from  that  of  the  ground  plane  and  that 
the  theory  and  measurements  are  in  close  agreement  Overlap 
or  hinge  surface-patch  dipoles  are  used  to  insure  continuity  of 
current  at  the  plate  to  plate  junction 

Fig.  7  shows  the  input  admittance  of  the  geometry  of  Fig  3 
with  the  addition  of  a  parasitic  wire  element  1  5  cm  above  the 
ground  plane  and  parallel  to  the  fed  monopole  The  admittance 
is  shown  as  a  function  of  the  separation  J  of  the  two  wires. 
Fig  8  shows  a  similar  geometry  except  that  the  parasitic 
element  now  touches  the  ground  plane 

Fig.  9  shows  the  measured  and  calculated  input  impedance 
of  a  T  bar  fed.  rectangular  cavity-hacked  slot  antenna  opening 
into  a  finite  ground  plane 

IV  CONCLUSION 

\  moment  method  solution  to  the  problem  of  radiation  or 
scattering  from  geometries  consisting  of  wires,  open  or  closed 
surtaxes,  and  wire  surface  junctions  has  heen  presented  The 
technique  is  a  Galerkin  solution  based  on  the  sinusoidal  reav 
tion  formulation  The  expansion  and  test  modes  were  chosen 
to  obtain  continuity  of  current  on  the  composite  structure, 
the  pr'ner  p/p  behavior  of  the  surface  current  density  at 
attachment  points,  and  facilitate  the  numerical  evaluation  of 
the  impedance  matrix  Numerical  examples  presented  illustrate 
the  accuracy,  versatility,  stability,  and  computational  efficiency 
of  the  method  Future  work  will  center  on  further  improving 
the  speed  and  computational  efficiency  of  the  method  and 
adding  the  presence  of  dielectrics 
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!  ig  7  Input  admittance  versus  separation  of  driven  mon  .pule  and 
parasitic  wue  element  0.0 1 S  m  above  the  0  0914  v  0  91 4  m  ground 
plane  Driven  monopole  height  =  0  431  m  Parasitic  wire  length 
0  842  m  /=  150  MHz 
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Fig  H  Input  admittance  versus  separation  of  driven  monopole  and 
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Near  fields  of  a  vector  electric  line  source  near  the  edge  of  a  wedge 
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Simple  closed  form  expressions  are  obtained  for  the  fields  of  an  electric  line  source  with  transverse 
or  longitudinal  current  flow  parallel  to  and  near  the  edge  of  a  perfectly  conducting  wedge  The 
effect  of  an  edge  which  is  not  perfectly  sharp  is  investigated  bv  considering  a  cylinder  tipped 
half  plane 


I  INTRODUCTION 

The  problem  of  scattering  by  wedges  has  been 
treated  extensively  by  many  authors  [Bowman  et 
at..  1969;  Felsen  and  Marcuvitz.  1973]  Usually  the 
analysis  involves  asymptotic  methods  for  evaluating 
far  fields,  with  little  emphasis  on  near-field  calcula¬ 
tions  It  is  the  near-field  aspect  with  which  the 
present  paper  is  concerned 

An  interest  in  the  near  zone  fields  of  a  source 
near  an  edge  arises  when  formulating  a  moment 
method  solution  to  the  problem  of  a  wire  antenna 
mounted  near  the  edge  of  a  plate  or  wedge.  In 
a  previous  paper  |,Vfnman  and  Pozar.  1978]  the 
authors  presented  a  moment  method  technique 
whereby  problems  with  wires  mounted  0.  I\  or 
greater  from  an  edge  may  be  treated.  The  analysis 
presented  here  forms  a  first  step  in  extending  this 
technique  by  synthesizing  an  expansion  mode  to 
treat  the  case  of  the  wire  near  an  edge  A  wire 
attached  near  the  edge  of  a  wedge  may  have  current 
transverse  or  parallel  to  the  edge  Thus  a  vector 
electric  line  source  parallel  to  the  edge  but  with 
current  transverse  or  parallel  to  the  edge  is  consid¬ 
ered  Simple  expressions  for  the  near  zone  fields 
are  derived  which  exhibit  not  only  the  correct  edge 
behavior  but  the  source  singularity  as  well 

2  TH KORY  H)R  PEREEC  TLY  CONDUCTING  WEDGE 

A  Introduction  The  geometry  of  the  problem 
is  shown  in  Figure  I  The  field  point  is  (p.  <b>  while 
the  electric  line  source  is  at  (p'.  <b'|.  0  £  <t>.  tb' 
£  2tt  —  <b„  The  line  source  is  parallel  to  the  edge 
of  the  wedge  and  has  components  /  ,  /  and  / 
of  current  in  the  p.  <i>.  and  z  directions,  respectively. 

<  .ipirtphi  ft  i hv  the  4 mem  ait  Geophysical  Union 


The  current  is  assumed  to  be  constant  with  respect 
to  z.  Thus  the  source  current  density  may  be  written 
as 

I 

J  (p,  <b)  =  (/„  p  +  /,<b  +  /.  ; )  —  6(p  p' ) 

P 

( I ) 

The  method  of  solution  is  to  use  the  two-dimen¬ 
sional  dyadic  Green's  function  for  the  wedge  to 
find  the  magnetic  field  due  to  the  vector  current 
J  Then,  since  the  field  and  source  points  are  close 
to  the  edge,  small  argument  approximations  are  used 
for  the  Bessel  functions  which  occur  in  the  resulting 
infinite  series.  The  series  is  then  summed  in  closed 
form,  preserving  the  source  singularity.  Maxwell's 
equations  are  used  to  find  the  E  fields  Of  course, 
with  the  fields  of  an  electric  source  known  the 
fields  of  a  magnetic  source  may  be  found  from 
duality. 

B  Two-dimensional  d\adic  Green's  Junction 
The  two-dimensional  dyadic  Green's  function  is  a 
solution  of  the  vector  wave  equation  [7a/.  1973] 

(V  x  V  *  -*;)G„.(R|R  )  =  V  x  [iSiR  -  R  i]  (2) 

where  k  =  2-rr/k,,.  I  is  the  unit  dyad.  R  =  pp 
<b<t>.  R  =  p'  p  +  <b  ’  <i>.  and  ,  satisfies  the  boundarv 
condition  n  x  V  x  G_,  =  6  on  S.  The  fields  are 
then  found  from 

HlR)  =  S  G„,  i R  R  T  -  J i ft  )  ds  ( 3 1 

and.  assuming  e'“"  time  dependence. 

I 

E(Rl  =  -  |T  x  HiRi  -  J(Ri|  (4) 

Juif  , 

The  integration  in  (3)  is  over  the  source  region 
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to  the  source  polarization  and  subscripts  refer  to 
the  field  component  Then,  from  ( 5 >-(7 ).  H  =  H  ‘ 
=  HZ  =  HZ  =  H\  =  0.  and 


CROSS-SECTION 


WEDGE  SURFACE  3 


H"  (p.  <t>l 


<  2  TT  - 


f  '  aKj,  <  A.p >7 ,  (Ap' ) 

i  (/a  > - cos  ud>  sin  rid 

Jo  rf,  u  +  Mv  -  *  > 


/ftp.  <i>)  - 


Hg  I  (iCumciry  of  ihe  line  source  near  an  edge 


Using  ( I )  in  (3)  gives 

HiRi  /„(!.  i  R 1 R ' )  *  p 

+■  /%G„  iR  R  )  •  ii>  f  /,  G„  i  R  ,R  I  •  i  (5 1 

The  Green's  function  G^.  is  given  by  Tai  [1971. 
1973]  i note  that  from  (I).  V  •  G„,  =  0.  so  (hat 
only  the  VI  and  N  vector  wave  functions  are  needed 
in  the  eigenfunction  expansion  of  G,„,): 

G„.  |R  K  > 

2  f  •  vt  v.  +  \  vi' 

\  dA  >  — .  ■— - (6) 

(2tt  d),.)  J(l  ,  (I  +  <V)(A'  -  k  ) 


(  2  TT  d>„) 

„  7,  t \p)  J\  tAp’  ) 

\  cos  ad)  cos  ud>  ld| 

( I  T  6„)(A"  -  k  ;) 


H  ip.  <t>) 


(2 tt  -  6„)p 


22,  a\ /,  (Ap)7,  (Ap’  I 

ii\  N - —  cos  od>  sin  ad> ’  (10) 

“  (I  +  o.,)< A'  -  k') 


H  .  Ip.  ■»)  = 


( 2  TT  —  <t>  ) 


JO  n  =  <» 


2,  \  '  J\[Kp)J  (Ap’l 

y  - t — — 

-7,  (I  +  —  k~) 


sin  a<t>  sm  ad)  (III 


f  I  n  =  0 

s»  =  y 

10  GO 


a  —  nv  —  n 


(2ir  -  d)„) 


n  =  0,  1.  2,  .  and  the  vector  wave  functions 

are  (the  two-dimensional  form  of  the  vector  wave 
functions  given  above  is  most  easily  found  from 
the  three-dimensional  form  given  by  [Tai.  1971] 
by  letting  h  ~  0) 

a  sin  cos 

M,  c  -  7t(Ap)  a<t>p  -  A/'.  <\p>  a«t><b 
i/*”  p  cos  sin 

(7) 

cos 

4  =  /.(Ap)  (i<t> : 

sm 

Consistent  with  Ten's  [  1 97 1 1  notation.  VI'  and 
V  in  (6)  refer  to  M  and  N  with  primed  spatial 
(source)  coordinates,  whereas  elsewhere  the  prime 
refers  to  differentiation  with  respect  to  the  argu¬ 
ment 

In  defining  field  components,  superscripts  refer 


The  integrations  in  ( 8 )— ( II)  can  be  evaluated  [Tai. 
1971]  to  give 

jirl  22,  a  cos  ad>  sin  ad) 

HU  p.  d>)  = - - -  >  - : - 

(2tt  -  d>„)p'  (i+8„) 

(  J,{kp)H;'{kp' )  p  <  p' 

{  ,  ,  (12) 
l  H  ;\kp)J,_  <* p' )  p  >  p 

Juki  2-,  cos  a<b  cos  ad)’ 

H'ip.  d»  = -  >  - 

< 2 -tt  -  d>„)  r~f,  (i  +  8„i 

f  J  i  kp )H  "  (A p ’ )  p  <  p 

(13) 

{  H  ,  (kp)J\  Ik  p' )  p  >  p 

- /  — /  _2,  a  ctss  a<t)  sir!  ad) 

H  |p.  'Si  -  \ 

1 2  TT  dv.lp  “  (  1  F  ,),) 


J,(kp)H  {k p  I  p<  p 


H  ,  Ik  p)7.  ( kp  I  P  p 


jixkl  —2,  sin  ad>  sin  ad)' 

H .  I  p.  d>)  =  - - S  — - 

(2  it  -  d>„)  „„  (I  Eft,,) 


f  (kp)H',  ‘(k  p') 
I  H  ,  ik  p  )J ,  (A  p 


)  P  ■  p 
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Up  tt)  this  point  the  field  expressions  are  completely 
rigorous 

<  Sear-field  calculation  The  source  and  field 
points  are  non  restricted  to  be  close  to  the  edge 
of  the  wedge,  so  k  p  <  I  and  kp  •'I  Using  small 
argument  approximations  for  the  Bessel  functions 
in  ( 12)  (15)  and  trigonometric  identities  sields 


H  ,,  <  p.  4) 


/  P 


:<Zn  4„> 

p  cos  i  (6  +•  4  )  p  p  cos  rl6  <t>  l  p 


//'  ( p.  <t> » 


<2tt  >5..  Ip' 


V  ” 


|sm  n  (6  v  4 ' )  sin  o(4  4' 


p  <  p 


1 25 1 


t  lr>) 


H " (p.  <i>) 


(2rr  4M)p‘ 


//„( p.  4) 

(2tt  -  4.  )p 


W„(p-  4) 


[sinutib  <b  )  sina'ib  4  )|  p  *s  p’ 

^5"  (  )  [cosa(4  +■  4' )  +  cosu(4  -  4’  )|  p  <  p 

1  V’  (*  V  . 

^  ^  —  J  (cost*  (<t>  -t-ft)  )  +  COS  (I  (<t>  <j)  )j  +•  I  p  >  p 

1  v/fV 

/ _  I  I  [sin a(4  -p  ib  1  sinnbb  —  ib  i| 

2  „  \p  / 

1  V  /p  V 

—  —  I  [sin  a  (4  +•  <b'  )  -  sina(4  4'  )| 

rS(5) 


p  <  p 


p  >  P 


1  1  H) 


[cos  a(4  p  4  )  -  cos  a  (4  4 


P  <  P 


(  2it  -  4  ,)p 


-!  A  /p'V 

- 2,  I  —  I  [cosa  (4  +  4' )  —  cos  a(4  - 

12  “  \  p  / 


( ls>) 


4' )[  p  >  p' 


Using  the  summation  formulas  (A3)  and  (A4) 
given  in  the  Appendix,  after  simplification,  the  where  v  -  -n/(2n  -  4,)  and 


H“(  p.  4)  - 


desired  closed  form  results 

f.p'p' 

2(2tt  4„) 

sin  i  (4  +•  4' )  sin  v(4  4  ) 

R.  R . 

/„?'■  1 

//  '  (p.  4)  =  — - 

2(2tt  4  ) 

p  -  p’  cos  1(4  +  4  )  p'  -  p  cos  o  4  -  4  ) 

R  R 


R .  =  p‘‘  +  p''  -  2p'  p'  ‘  cos  i<4  i  4'  I 

The  E  fields  could  be  found  either  by  using  (4) 
on  ( 1 2 )  ( 15)  and  applying  the  procedure  of  small 
argument  approximations  and  summation  of  the 
(20)  senes  or  by  using  (4)  on  (20h23)  directly.  Identical 
results  are  obtained  with  either  method  for  all  E 
components  except  for  E  .  which  would  reduce 
to  zero  if  the  curl  of  (2-2).  (23)  were  taken  Thus 
I  is  found  from  i4|  and  ( 14)  1 15).  while  the  other 
E  components  are  found  by  using  (4)  and  (20).  (21 ) 
After  simplification.  £*  =  £"  =  £  =  £  =  ij. 

and 


H  t  p.  'b) 


/  p  1  p 

2(2—  4  1 


sin  i  (4  f  4’ )  sm  i(4  -  4' ) 

R  R 


(21) 


(22) 


£''(p.  4)  = 


/„  Wpp' ) 

2/<ut  , ( 2 rr  4  ,) 

p  p  (p  +p  )  cos  i  (4  +  4  1 
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From  (4),  (14).  and  ( 15). 
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Using  small  argument  approximations  in  the  Bes¬ 
sel  f  unctions  in  (28)  gives 
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Using  the  summation  formula  i  A7)  of  the  Appendix 
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Equations  ( 20)— ( 23 ).  ( 24)— < 27).  and  (30)  are  the 
complete  Held  expressions  for  an  arbitrary  vector 
electric  line  source  and  are  rigorous  in  the  near-field 
limit,  including  the  proper  source  singularity  Note 
the  proper  edge  behavior  in  the  fields  as  p  -•  fl 
Also  note  that  the  expressions  for  H  in  (20)— (23 ) 
are  independent  of  frequency,  indicating  a  quasi- 
static  result  Hence  (20)- (23)  could  have  been  ob¬ 
tained  from  a  conformal  mapping  approach 

Figure  2  shows  the  magnitude  of  E*  versus 
distance  along  the  surface  of  the  wedge  for  various 
wedge  angles  and  k  -  I  m.  p'  =  0 .1  m.  A'  = 
0.  and  =  I 

Although  the  transverse  Field  components  are 
singular  at  the  edge  of  the  ideal  wedge,  no  such 
singularity  exists  on  any  physically  realizable  elec¬ 
tromagnetic  scatterer  The  reason  for  this  difference 
is  that  the  ideal  wedge  has  a  perfectly  sharp  edge 
with  a  discontinuous  normal  vector,  u  hile  in  phx  si 
cal  realitv  any  real  edge'  will  ha\e  a  small  radius 
of  curvature  so  that  the  normal  vector  at  the  edge' 
will  change  ..ontnuouxly  as  it  moves  around  the 
'edge  '  This  fact  w3s  recently  used  by  Rhodes 
(1471 1  in  connection  with  physically  realizable 
antenna  aperture  distributions 
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Pig  2  Magnitude  of  E*  (phase  Xf)  versus  normalized 
distanve  along  huih  sides  of  the  wedge  surface  for  various  wedge 
angles  p  0  I  m  *  I  m  *'  0.  /,  I 

Another  interesting  observation  is  that  some  field 
components  away  from  the  edge  become  singular 
when  the  line  source  approaches  the  edge  (p'  — » 
0)  This  is  ir,  agreement  with  reciprocity.  Of  course, 
this  situation  could  never  happen  in  practice  because 
the  perfectly  sharp  wedge  does  not  exist 

In  order  to  get  a  feeling  for  the  field  behavior 
at  a  physical  'edge'  the  mathematically  ideal  edge 
of  the  wedge  of  this  section  is  removed  by  adding 
a  circular  cylinder  tip  Although  this  geometry  is 
not  an  exact  model  of  a  physical  edge,  it  should 
allow  a  qualitative  view  of  the  field  behavior,  with 
the  advantage  that  it  is  mathematically  tractable 

3  THE  CYLINDER  TIPPED  HALF  PLANE 

By  using  the  same  method  as  that  outlined  in 
section  2,  the  closed  form  field  equations  for  a 
cylinder  tipped  wedge  may  be  obtained  Here,  only 
the  expression  for  F *  for  the  special  case  of  a 
half  plane  (6,  =  Of  with  an  /,  line  source  1/ 

/,  =  0)  on  the  surface  of  the  half  plane  (Pi’ 

0)  is  presented 

The  geometry  for  the  cvlinder  tipped  hall  plane 
is  shown  in  I  noire  '  (he  two-dimensional  dvudiv 
Green  s  function  mav  be  obtained  trom  the  three- 
dimensional  result  given  by  Tut  ||97||  as  described 
in  section  2  The  Greens  function  is  seen  to  consist 
of  a  sum  of  the  half-plane  Green's  function  and 
terms  which  account  tor  scattering  from  the  cy¬ 
lindrical  tip  Thus 


I  ig  5  Geomelry  of  an  electric  line  source  on  the  surface 
and  close  lo  the  edge  of  a  cshnder-ttpped  half  plane  Current 
polarization  is  in  the  <t»  *»r  1 1  direction 
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w  here  F.*'  *r  ‘1‘"  ( R  )  is  given  by  (26)  and  (27)  Omit¬ 
ting  the  detailed  derivation,  for  k p  <  I  and  k p 
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(31) 

The  total  F\  field  for  t  =  0.  6  -  0.  6.,  =  Q. 
p'  =  0  I  m,  \  =  I  m  is  plotted  versus  p  =  x  for 
v  artous  values  of  a  for  the  cy  linder-ttpped  half  plane 
in  Figure  4.  Also  shown  in  this  figure  is  the  F 
field  for  a  perfectly  sharp  half  plane,  which  would 
correspond  to  a  0  (no  cylinder  tip >  Figure  5 
shows  the  same  curves  for  p'  =  0.01  (A  =  I  mi 
The  p  -  t  axis  is  logarithmic  to  give  an  expanded 
view  of  the  edge  area  First  of  all.  note  that  for 
a  =  0  (no  cylinder)  the  field  is  becoming  singular 
as  expected  Now.  for  a  small  but  finite  cylinder 
radius  a  the  F  field  increases  as  the  edge  is 
approached  but  reaches  a  maximum  and  goes  to 
zero  at  p  a.  as  is  required  bv  bound.irv  conditions 
It  can  he  seen  trom  Figures  4  and  5  that  :he  held 
reaches  a  higher  maximum  for  smaller  j  but  is  not 
singular  except  for  a  -  0. 

Figure  6  shows  the  total  F  field  for  i  0. 
4>  =  0.  6,,  =  0.  p'  =  0  I  m  \  I  m.  and  a 

5  <  10  1  m  versus  p  =  x  on  a  rectangular  scale 
for  the  ideal  half  plane  and  the  cylinder-tipped  half 
plane,  along  with  a  sketch  of  the  relative  size  ot 
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f  tg  4  Magnitude  of  t.\  (phase  9(f)  versus  position  (p 
x )  along  the  surface  of  a  cylinder-tipped  half  plane  for  various 
values  of  cylinder  tip  radius  a.  p'  =  0  I  m.  x  =  1  m.  d>' 

0.  <t>  0,  I.  6,  0  (half  plane) 


f  ig  6  Magnitude  of  E\  (phase  90°)  versus  position  (p 
t)  along  the  surface  of  a  cylinder-tipped  half  plane  compared 
with  that  for  an  ideal  half  plane  The  relative  si/.c  of  the  cylindrical 
tip  is  also  shown,  p  =  0.1  m,  k  -  1  m,  6'  =  0.  6  =  0. 
/_  =  I .  a  =  p'/2(J0 


Fig  5  Magnitude  of  f  \  (phase  'Xf)  versus  position  (p 
x)  along  the  surtave  of  j  cylinder-tipped  half  plane  for  various 
values  of  cylinder  tip  radius  a.  p'  =  0  01  m.  X  =  I  m. 

-  0.  4>  0.  I.  0  (half  plane) 


the  cylinder.  Figure  7  shows  the  same  curves  for 
a  =  5  x  10  m 

4  CONCLUSION 

Simple  closed  form  expressions  have  been 
derived  and  presented  for  the  fields  of  an  electric 
line  source  with  arbitrary  current  flow  near  and 
parallel  to  the  edge  of  a  perfectly  conducting  w  edge 
The  effect  of  a  physical  edge  (not  perfectly  sharp) 
has  been  investigated  by  comparing  the  lields  near 
the  end  of  a  cylinder-upped  half  plane  with  those 
near  ihe  edge  of  an  ideal  half  plane 

The  above  results  show  that  for  even  a  small 
departure  (torn  an  ideal  edge  the  fields  will  not 
be  singular  Since  anv  physical  scattering  edge  has 
a  certain  degree  ol  roundness  or  roughness,  singu¬ 
larities  will  not  be  found  in  the  edge  vicinity  It 
should  be  mentioned  that  there  is  not  universal 
agreement  on  the  interpretation  of  the  dilterences 
between  ideal  edges  and  phy  sical  edges  [l  ee.  I ^721 
However,  it  is  clear  from  the  data  presented  that 
for  a  source  in  the  near  vicinity  of  an  edge  the 
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Z  3  *  10' 


E  quating  real  and  imaginary  parts  gives  the  desired 
results 


Fig  7  Magnitude  of  f.\  iphase  =  -9(f)  versus  position  (p 
c)  along  the  surface  of  a  cylinder-tipped  half  plane  compared 
with  that  for  an  ideal  had  plane  The  relative  size  of  the  cylindrical 
up  is  also  shown,  p  =  0  I  m,  \  =  I  m.  6'  =  0.  <6  =  0. 
/.  =  l.o  p'/2000 


‘edge’  behavior  of  the  fields  is  confirmed  to  a  region 
a  few  thousandths  of  t  wavelength  from  the  edge 


APPENDIX  SUMMATION  OF  SERIES 

Starting  with  the  known  result 
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Replacing  a  by  ae*  and  equating  real  parts  gives 
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formulation  are  speed,  accuracy,  convergence,  and  versatility.  Techniques 
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based  on  the  piecewise  sinusoidal  reaction  formulation. 
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I.  INTRODUCTION 

In  a  previous  paper  LlJ  the  authors  presented  a  moment  method 
solution  for  wires,  rectangular  plates,  and  wire/plate  attachments. 

The  purpose  of  this  paper  is  to  present  some  details  of  the  formulation 
and  of  the  computer  code,  and  especially  to  show  how  the  choice  of  in¬ 
tegral  equation  and  expansion  and  testing  functions  impact  on  the  ver¬ 
satility,  accuracy,  computational  efficiency,  and  ease  of  use  of  this 
and  other  similar  formulations. 

Briefly,  the  solution  [1]  is  based  on  the  reaction  integral  equation 
and  employs  speed-of- 1 ight  piecewise-sinusoidal  (PWS)  wire,  surface, 
and  attachment  dipole  expansion  modes  which  are  placed  on  the  composite 
w>e/surface  geometry  in  an  overlapping  array  such  that  continuity  of 
current  is  enforced.  The  weighting  or  test  modes  are  chosen  identical 
to  the  expansion  modes,  and  thus  the  solution  is  a  Galerkin  method, 
yielding  a  symmetric  impedance  matrix.  Electric  test  sources  are  used, 
so  open  as  well  as  closed  surfaces  can  be  treated. 

II.  MODE  LAYOUT 

The  mode  layout  will  be  described  with  the  aid  of  an  example, 
a  monopole  on  a  bent  plate  as  shown  in  Figure  1.  The  details  of  the 
modes  have  been  given  previously  [l  ].  Since  this  problem  involves  wire, 
plate,  and  attachment  modes,  the  impedance  matrix, [  Z ],  can  be  symbol¬ 
ically  shown  as: 
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w/w 
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*/p 

"A 

p/w 

■ 

p/p 

% 

A/w 

A/p 

W  =  WIRE 
P  =  PLATE 
A  =  ATTACHMENT 


W i re  Modes 

Figure  ?a  shows  the  wire  broken  into  three  segments,  or  two  PWS 
dipole  modes  L?J.  Note  that  while  the  wire  current  generally  does  not 
vanish  at  the  attachment  point,  the  wire  mode  current  does.  Thus  the 
need  for  an  attachment  mode,  described  below. 

Surface  Modes 

in  the  computer  code  a  rectangular  plate  is  defined  by  specifying 
the  coordinates  of  three  consecutive  corners,  and  the  plate  segmentation 
in  the  two  orthogonal  directions.  Figure  2b  shows  plate  1  of  the  bent 
plate  geometry  divided  into  three  segments  in  the  l-to-2  direction  and 
three  segments  in  the  2-to-3  direction.  Six  overlapping  PWS  surface 
dipoles  are  shown  as  arrows  in  the  l-to-2  direction,  and  six  in  the 
?-Lo-3  direction.  Not  shown  are  the  modes  on  plate  2. 

Attachment  Modes 

The  purpose  of  the  attachment  mode  is  to  establish  continuity 
of  current  at  the  wire/plate  junction  and  also  to  insure  the  proper 
I h  singularity  of  the  plate  surface  current  density  in  the  vicinity 
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of  the  attachment  point.  Figure  2c  shows  the  PWS  attachment  dipole 
consisting  of  a  circular  disk  monopole  with  1/:  surface  current  density, 
and  a  wire  monopole  with  non-zero  current  at  the  attachment  point. 

The  attachment  point  must  be  at  least  0.1>  from  the  edge  of  a  plate, 

but  it  need  not  be  at  the  center,  at  a  corner,  or  have  any  special  location 

with  respect  to  the  surface  dipole  modes. 

Overlap  Modes 

When  two  plates  intersect  additional  surf  ace-patch  dipole  modes 
are  required  to  allow  a  nonzero  continuous  current  at  the  plate-to-pl ate 
junction.  These  modes  are  termed  overlap  modes.  They  are  identical 
to  the  surface  dipoles  on  plates  1  or  2,  except  that  the  dipole  angle 
may  differ  from  180°.  Figure  2d  shows  two  overlap  modes  at  the  junction 
of  the  bent  plate.  The  edges  of  the  overlap  surface  dipoles  need  not 
coincide  with  the  edges  of  the  surface  dipole  modes  on  plate  1  or  2, 
thus  allowing  the  intersection  of  plates  of  different  si2e.  The  code 
automatically  checks  to  see  which  plates  intersect  and  inserts  the  over¬ 
lap  modes. 

Toep I i tz-L ike  Properties 

The  plate-to-plate  block  of  the  impedance  matrix  is  shown  below: 


p'/p, 

Vp, 

Vo 

Vp, 

Vp2 

Vo 

°/p, 

°/p. 

% 

P,  =  PLATE  I 
P2  =  PLATE  2 

0  =  OVERLAP 
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The  various  partitions  show  plate  1-to-platr  1  surf  .'cc  dipole  mutual 
impedances,  plate  1-to-plat.e  ?  mutuals,  overlap-to-plate  1  mutuals, 
etc.  Hue  to  the  regular  nature  of  the  placement  of  surface  dipole  modes 
on  a  plate,  considerable  savings  in  computation  time  is  available  in 
computing  the  plate  1-to-plate  1  block  and  the  plate  2-to-plate  2  block. 

For  example,  consider  the  plate  1-to-plate  1  block,  for  which  the  modal 
layout  is  shown  in  Figure  2b.  Note  that  Z,  0  =  Z,  .,  Z.  .  =  Z_ 

i  )  u  j  T  JjD 

Zj  ^  =  -  Zq  g,  etc.  On  an  arbitrary  plate  it  is  only  necessary  to 
compute  the  mutual  impedances  between  the  first  mode  in  the  l-to-2  direction 
and  all  the  modes  on  the  plate,  and  the  mutuals  between  the  first  mode 
in  the  2-to-3  direction  and  all  2-to-3  modes  on  the  plate.  The  detail 
of  the  plate  1-to-plate  1  block  is  shown  below: 


Fiie  x  represents  elements  which  must  be  computed,  and  the  0  represents 
elements  which  can  be  obtained  from  the  toepl itz-1 ike  properties.  If 
on  a  given  plate  there  are  modes  in  the  l-to-2  direction  and  N?3 
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modes  in  the  P-3  direction,  then  only  N ^  +  ?  N^  of  the  (N^  +  N^)" 

elements  need  to  be  computed. 

III.  COMPUTATIONAL  ADVANTAGES  OF  MODE  CHOICE 
Convergence 

A  rapidly  convergent  set  of  modes  will  minimize  the  size  and  thus 
the  time  reguired  to  compute  the  impedance  matrix  and  the  required  storage. 
Figure  3  shows  a  comparison  of  the  input  admittance  of  a  monopole  on 
a  small  plate  computed  using  15  PWS  modes  to  that  using  68  pulse  modes  [3]. 
Noting  that  15  PWS  modes  show  better  agreement  with  measurements  than  68 
pulse  modes,  it  is  clear  that  the  PWS  basis  converges  relatively  fast. 
Computation  of  Impedance  Elements 

Two  algebraically  equivalent,  but  numerically  different,  techniques 
for  evaluating  the  elements  will  be  presented,  and  their  relative  ad¬ 
vantages  in  computing  various  elements  described.  The  discussion  below 
will  concentrate  on  surface-to-surface  dipole  and  di sk-to-surf ace  dipole 
impedances  which  occupy  the  majority  of  the  time  in  the  computation 
of  the  upper  triangular  part  of  the  symmetric  impedance  matrix. 

A  general  impedance  matrix  element  Z  is  defined  as  [ l] 

mn 


where  e^e^  are  independent  coordinates  on  the  surface  of  the  n-th  ex¬ 
pansion  mode.  Jn  is  the  current  density  of  the  n-th  expansion  mode, 
■and  Fm  is  the  free-space  electric  field  of  the  m-th  test  mode,  which 
can  he  written  as 
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(2) 


£m^el,e2^  (  J  S0(e1,e2;t1,t2).  JrT)(t1,t2)dt1dt2, 

]  l2 

whore  J|n  is  the  current  density  of  the  m-th  test  mode,  t-j  and  t2  are 

independent  coordinates  on  the  surface  of  the  m-th  test  mode,  and  §Q 

is  the  free-space  Green's  function.  Equations  (1)  and  (2)  show  that 

an  impedance  element  in  general  requires  a  4-fold  integration,  i.e., 

two  to  find  Em  and  two  to  integrate  over  the  surface  of  the  n-th  expansion 

inode.  A  useful  property  of  the  PWS  modes  is  that  E  is  known  in  closed 

m 

form  for  a  wire  [4]  or  a  surface  [5.1  monopole,  and  requires  only  one 
simple  numerical  integration  for  a  disk  monopole.  Each  dipole-to-dipole 
Z  is  the  sum  of  four  monopole-to-monopole  impedances. 

The  second  way  of  computing  the  impedance  elements  is  to  consider 
the  surface-dipole  mode  current  and  disk  monopole  current  as  being  made 
up  of  PWS  filaments.  The  order  of  integration  in  Equations  (1)  and  (2) 
can  be  interchanged  to  give 


Zmn  =  /  i  We2’fc2)de2dt2’ 
e2  l2 

where  *inn(e2,t2)  represents  the  mutual  impedance  between  two  filaments 
lying  in  e1  and  t^  directions.  That  is. 


z  (e 
mnv 


‘Vel’e2^  r  'Vei»e2;t1,t2).Jm(t1,t2)dt1de1 

e ,  t , 


(4) 


An  advantage  of  the  PWS  modes  is  that  the  zmn(e2»t2)  of  Equation  (4) 
are  known  in  closed  form  ^6i. 
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* 


Ar:  v1v.int.aqo  of  the  second  method,  similar  to  that  employed  by 
Popovir  and  PopovicL?!  occurs  when  computing  the  impedances  between 
two  surface  patch  monopoles  if  either  the  surface  current  directions 
are  parallel,  or  the  vector  transverse  to  the  surface  monopoles  are 
parallel.  If  the  expansion  and  test  surface  monopoles  are  each  repre- 

O 

son  tod  by  M  filaments,  then  Equation  (3)  requires  M  evaluations  of 

/  (of).  However,  no  more  than  PM  of  the  z  are  different.  By 

inn  /  /  mn 

storing  the  PM  values  prior  to  computing  the  mutual  impedance,  the 
computation  time  will  be  proportional  to  2M  rather  than  M  . 

Another  time-saving  technique  deals  with  the  slow  convergence 
of  the  numerical  integration  of  Equation  (3)  when  computing  the  mutual 
impedance  of  two  touching  surface  monopoles.  This  slow  convergence 
is  a  result  of  the  fact  that  the  imaginary  part  of  the  mutual  impedance 
of  two  PWS  filaments  has  a  logarithmic  singularity  as  the  distance  between 
the  two  filaments  gets  small  [8l  For  small  separation  x,  the  reactance 
between  two  filaments  can  be  written  as  L 8  ] 

X(x)  =  Cj  +  C?  ln( x) .  (5) 

The  constants  and  can  be  easily  evaluated  and  the  logarithmic 
singularity  integrated  analytical ly.  Figure  4  shows  the  self-reactance 
of  a  surface  monopole  computed  using  Equation  (1),  Equation  (3)  directly, 
and  fgii.it  ion  (3)  with  the  above  two  techniques  (i.e.,  only  M  evaluations 
of  /  and  extracting  logarithmic  singularity).  Also  shown  is  the  com¬ 
putation  time  for  each  of  the  three  methods;  all  times  shown  in  this 
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r  .ire  lor  the  Datacraf t  60?4  which  is  about  an  order  of  magnitude 
slower  than  an  IBM370-165.  In  using  Equations  (1)  and  (3)  directly 
the  surface  monopoles  are  separated  by  10"4*  to  avoid  the  singularity. 

No  separation  is  required  if  the  logarithmic  singularity  is  removed. 

From  Figure  4  it  can  be  seen  that  not  only  is  the  use  of  Equation  (3) 
with  the  above  discussed  two  techniques  faster  than  the  other  methods, 
but  it  is  also  more  accurate. 

One  case  where  Equation  (1)  has  an  advantage  occurs  when  computing 
the  mutual  impedance  between  a  disk  monopole  (of  an  attachment  mode) 
and  the  set  of  surface  dipoles  on  a  plate  which  is  parallel  to  the  plane 
of  the  disk.  Here,  advantage  can  be  taken  of  the  fact  that  only  the 
E  field  of  the  disk  monopole  is  needed  (since  E.=0  and  E  z.J  =0), 
and  also  E  is  independent  of  <$>.  Thus,  much  time  can  be  saved  by  filling 
a  table  of  values  of  E^  versus  p.  This  table  can  then  be  interpolated 
to  find  Em  on  every  surface  dipole  mode  on  the  plate,  making  the  evalu¬ 
ation  of  Equation  (1)  very  fast.  An  order  of  magnitude  saving  in  com¬ 
puter  time  is  then  realized  for  small  plates,  with  substantially  greater 
savings  for  larger  plates. 

V.  CONCLUSION 

This  paper  has  demonstrated  that  in  a  moment  method  surface/wire 
formulation,  the  piecewise  sinusoidal  modes  have  significant  advantages 
from  the  standpoints  of  convergence  and  the  fast  and  accurate  computation 
of  the  impedance  matrix.  Also,  techniques  for  reducing  the  run  time 
of  surface/wire  codes  were  presented.  To  show  how  the  speed  of  a  typical 
solution  can  be  improved  by  these  methods,  the  problem  of  a  \/4  monopole 
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mounted  on  a  l.Ox  square  plate  was  run  without  using  any  of  l he  techniques 
mentioned  in  this  paper,  and  again  using  all  of  the  mentioned  tTjiu  iques . 
Using  the  techniques  uecreased  the  run  time  by  a  factor  of  twenty  to  115 
sec.  on  the  Datacraft  6Q?4  (about  20  sec.  on  IBM  370-'lo5)>  with  an  incr- 
in  accuracy  (both  runs  used  the  symmetry  of  the  impede  •  matrix). 

Acknowl edgment 

The  authors  would  like  to  thank  Professor  J.H.  Richmond  for  supplying 
invaluable  suggestions,  advice  and  computer  subroutines.  Helpful  dis¬ 
cussions  with  Dr.  W.D.  Burnside  and  Dr.  N.W.  Wang  are  also  appreciated. 


35 


REFERENCES 


I.  E.  H.  Newman  and  D.  M.  Pozar,  "E lectromagnet i c  Modeling  of  Composite 
Wire  and  Surface  Geometries,"  IEEE  AP-26,  pp.  784-739,  Novemoer 
1979, 

y.  Y.  T.  Lin  and  J.  H.  Richmond,  "EM  Modeling  of  Aircraft  at  Low  Fre¬ 
quencies,"  IEEE  AP-23,  pp.  53-56,  January  1975. 

3.  A.  W.  Glisson,  "On  The  Development  of  Numerical  Techniques  for 
Treating  Arb i trari ly-shaped  Surfaces,"  Dissertation,  University 
of  Mississippi,  June  1973. 

4.  S.  A.  Schelkunoff  and  H.  T.  Friis,  Antennas:  Theory  and  Practice, 

Wi  ley,  1952. 

5.  J.  H.  Richmond,  D.  M.  Pozar  and  E.  H.  Newman,  "Rigorous  Near-zone 
Field  Expressions  for  Rectangular  Sinusoidal  Monopole,"  IEEE  AP- 
26,  p.  509,  May  1973. 

6.  J.  H.  Richmond  and  N.  F).  Geary,  "Mutual  Impedance  of  Non-planar- 
skew  Sinusoidal  Dipoles,"  IEEE  AP-23,  May  1975. 

7.  Z.  D.  Popovic  and  B.  0.  Popovic,  "Transformation  of  Double  Integrals 
Appearing  in  Variational  Formulation  of  Cylindrical  Antenna  Problems, 
publications  De  la  Faculte  Dielectrotechnique  de  L'Universite 

A.,  Belgrade,  Serie:  Electronique,  Telecommunications  Automatique, 

No.  64,  1971. 

■9.  W.  A.  Imbriale  and  P.  G.  Ingerson,  "On  Numerical  Convergence  of 

Moment  Solutions  of  Moderately  Thick  Wire  Antennas  Using  Sinusoidal 
Basis  Functions,"  IEEE  AP-21,  pp.  363-366,  May  1973. 


36 


(d) 


1  iqure  ?.  Layout  of  the  (a)  wire,  (b)  surface  dipole, 
(c)  attachment,  and  (d)  overlap  modes  for  a 
monopole  on  a  bent  plate. 


SELF-REACTANCE  (OHMS) 


Figure  A.  Self-reactance  and  computation  time  For  a 
.15X  by  .2X  surface  monopole  computed 
three  ways. 
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ANALYSIS  OF  A  MONOPOLE  MOUNTED  NEAR  OR  AT 
THE  EDGE  OF  A  HALF-PLANE* 
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ABSTRACT 

The  problem  of  a  monopole  antenna  mounted  near  or  at  the  edge  of  a 
half-plane  is  solved  using  a  hybrid  moment  method  solution  with  an  integral 
expression  for  the  exact  half-plane  Green's  function.  Results  are  presented 
for  input  impedance  of  a  X/4  monopole  vs.  distance  from  the  half-plane  edge, 
and  vs.  angle  for  an  edge-mounted  antenna.  Simple  expressions  for  the  surface 
current  density  induced  on  the  halfplane  are  derived  and  used  to  synthesize 
an  attachment  mode  for  use  with  the  author's  previously  developed  surface 
patch  moment  method  solution  when  a  wire  is  attached  at  or  near  an  edge  of 
a  finite  plate. 
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DAA629-76-G-0331  between  the  U.  S.  Army  Research  Office,  Research 
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Foundation,  Columbus,  Ohio 
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INTRODUCTION 


A  problem  of  continued  interest  is  that  of  antennas  mounted  on  finite 
conducting  bodies.  Antennas  on  spherical  bodies  have  been  analyzed^,  as 
*<■'1  I  as  antennas  on  more  general  structures^ ’ J .  The  authors  have  previously* 
presented  a  surface  patch  moment  method  solution  for  geometries  involving 
wires  and  rectangular  plates,  including  plate/plate  and  wire/plate  junctions, 
where  the  wire/plate  attachment  point  was  restricted  to  be  . 1 \  or  more  from 
a  plate  edge.  This  paper  presents  a  new  wire/plate  attachment  mode  that 
removes  this  restriction,  allowing  a  wire  to  be  attached  close  to  or  on  a 
half-plane  type  edge.  The  development  of  this  new  surface  patch  moment  method 
expansion  mode  is  based  on  the  canonical  solution  of  a  vertical  monopole 
mounted  near  or  at  the  edge  of  a  perf ectly-conduct ing  half-plane.  T!vs  canonical 
solution  is  a  hybrid  moment  method  using  the  exact  half-plane  Green's  function 
in  integral  form,  with  the  source  and  image  terms  separated  from  the  dif¬ 
fracted  field.  This  form  of  the  Green's  function  also  yields  the  near  zone 
surface  current  density  on  the  half-plane  in  the  vicinity  of  the  edge  and 
source,  which  is  used  to  synthesize  the  above-mentioned  surface  patch  ex¬ 
pansion  mode  for  attachments  near  or  at  an  edge.  °resented  are  graphs  of 
input  impedances  for  a  vertical  X/4  monopole  vs.  distance  from  a  half-plane 
edge,  and  for  a  \/4  monopole  mounted  on  a  half-plane  edge  vs.  azimuth  angle. 
Far-field  results  for  antennas  near  half-plane  edges  are  available  in'. 

The  problem  of  a  wire  antenna  in  free  space  can  be  solved  via  the  moment 
method  solution  of  an  integral  equation  using  the  free  space  Green’s  function^. 

If  the  wire  antenna  is  in  the  presence  of  another  scattering  structure  there 
are  two  possible  modifications.  The  first  is  to  replace  the  free-space  Green’s 
function  by  the  special  Green’s  function  for  the  particular  scattering  structure. 
The  second  approach  is  to  retain  the  free-space  Green's  function,  but  suitably 
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expand  the  currents  on  the  wire  and  the  scattering  structure,  and  solve  for 
the  expansion  coefficients  in  the  MM  solution.  The  advantage  of  toe  first 
approach  is  that  no  new  unknowns  are  introduced  into  the  MM  solution  and 
the  si/e  of  the  impedance  matrix  does  not  increase  over  that  of  the  isolated 
wire,  while  the  disadvantage  is  that  the  Green's  function  of  the  scatterer 
must  be  known  in  a  computationally  efficient  form.  The  advantage  of  the 
second  approach  is  that  the  Green's  function  of  the  scatterer  need  not  be 
known,  while  the  disadvantage  is  that  new  unknowns  are  introduced  and  the 
si/e  of  the  required  impedance  matrix  can  increase  substantial ly  if  the  scat- 

9 

terer  is  electrically  large.  (Note  that  the  MM/GTD  hybrid  method  can  be 
considered  to  be  of  the  first  type,  where  the  exact  Green's  function  is  replaced 
by  the  GTD  or  high-frequency  approximation.)  Both  approaches  are  used  in 
this  paper. 

II.  SOLUTION  FOR  VERTICAL  MONOPOLE  NEAR  EDGE  OF 
HALF-PLANE  USING  EXACT  GREEN'S  FUNCTION 


Theory 

The  geometry  of  the  problem  is  shown  in  Figure  1.  The  wire  and  half¬ 
plane  are  considered  to  be  perfectly-conducting.  The  basic  idea  of  the 
solution  is  to  expand  the  wire  currents  using  piecewise-sinusoidal  modes, 
find  the  fields  of  test  modes  (also  PWS)  in  the  presence  of  the  half-plane, 
and  compute  an  impedance  matrix: 


Zij=-J  Ij(y>  /  Ii(y0)Ey  (y|y0)dy0dy  (1) 

yj  yoi 

where  I^y)  is  the  PWS  test/expansion  mode  current, 

E  (y|y0)  is  the  half-plane  Greek's  function,  representing 
y  the  E  field  due  to  a  y  infinitesimal  electric 

dipole  of  unit  moment. 


43 


Since  the  test  and  expansion  modes  are  equivalent  we  have  a  Galerxin 
solution,  and  Zjj-Zj-.  The  wire  currents  and  input  impedance  then  follow 
from  the  solution  of  a  matrix  equation  of  the  form  zj  f ij  =  fvj  ,  where  [zj 
is  the  impedance  matrix,  [ I  is  the  column  vector  of  current  expansion  coef¬ 
ficients,  and  : vj  is  the  vector  of  mode  voltages. 

Thus,  this  solution  depends  upon  the  availability  of  the  half-plane 
Green's  function  in  a  form  suitable  for  fast  and  accurate  evaluation.  It 
was  first  attempted  to  use  the  familiar  form  of  the  half-plane  Green's  function 
expressed  as  an  eigenfunction  expansion,  such  as  given  by  T a i ^ .  This  repre¬ 
sentation  was  found  to  be  unsuitable  for  two  reasons.  First,  the  eigenfunction 
expansion  contains  an  infinite  series  of  infinite  integrals  that  contain 
source  and  image  s ingul ari ties  that  would  make  evaluation  of  self  or  over¬ 
lapping  mutual  impedance  terms  difficult  to  do  accurately.  Second,  the  eigen¬ 
function  representat ion  does  not  lend  itself  to  small  argument  approximations 
for  evaulation  of  near  zone  surface  currents  in  simple  closed  form.  It  was 
found  that  the  first  term  of  the  series  would  give  the  correct  edge  behavior, 
but  not  the  correct  source  singularity.  For  these  reasons,  an  integral  ex¬ 
pression  for  the  half-plane  Green's  function  given  by  Senior^  was  used. 

7  ^ 

From',  the  total  field  due  to  a  y  directed  infinitesimal  electric 
dipole  of  moment  \l  in  the  presence  of  a  conducting  half-plane  is  (e'l'jt  time 
dependence) 


(2) 


where 

rn 

-  /  H|'  ^ ( kR  cosh  u ) du  (3) 
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Defining  the  modified  Hertz  potentials 


dR 


j,  f  (kR  cosh  u)du; 


(12) 


mR 


TdR  =  ^  /  H^(kR'  cosh  u)du, 

Mr. 


(13) 


and  using  the  relation 


7 


J  H^(kR  cosh  u)du  =  2j  — 


■  jkR 


kR 


(14) 


and  the  fact  that  and  4^-  =  4^-  ,  allows  the  diffracted  field 

q  q  ^y 

to  be  written  as 


rdif  .  ^0 
■y  4ttj 


^  f>\  f  32 


9'k  cos  |  cos  -j  H^^(kR-|)  - 
/PPn 


3yayr 


k2  „dR 


*f-£—  +  k2l 


\  3yayr 


(15) 


where  E  =  Einc  f  Eref  +  Edif  =  E9'0’  ♦  Edif 

y  y  y  y  y  y 

dR  dR ' 

Note  that,  using  Equation  (14),  it  and  it  can  be  expressed  in  terms  of 
finite  integrals: 


"dR  *  - \  "  i  f*  HS?,(kR  cosh  u)du, 

o 


(16) 


dR-  _  i  e-jkR'  j  y  (2) 


?  e_kft ■ - ^  /  Hp'(kR'  cosh  u)du. 


(17) 
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■'  ■  ■■  i—h:"  mu  l,i  ii  ■  .ijijBn, 


Defining  the  modified  Hertz  potentials 


idR  :  ^  f  H^(kR  cosh  u)du; 
^R 


(12) 


dR'  l 


=  {  f  (kR *  cosh  u)du, 

Mr> 


(13) 


and  using  the  relation7 


'O 

/  Hp^(kR  cosh  u)du  =  2j  — 


-  jkR 


kR 


(14) 


and  the  fact  that  and  -1?~  =  4?—  ,  allows  the  diffracted  field 

dyQ  °y  oyQ  ay 

to  be  written  as 


-d  if 


4  TT  j 


<t> 


^  cos  cos  —2  ( k R  i )  - 


-jk 

— COS  -k  COS  —X  H* 

, -  d  do 

rpp, 


3y3yr 


.  2)  dR 

-  k  )  TT 


/  a2 
+  — a-- 


^yay 


-+k2  ,dR' 


(15) 


where  E  =  Einc  *  Eref  +  Edif  =  E9’0*  *  rdif 

y  y  y  y  y  ^y  ' 

dR  dR 1 

Note  that,  using  Equation  (14),  tt  and  tt  can  be  expressed  in  terms  of 
finite  integrals: 


dR  1  e~jkR  j  ?  u( 2) , . 0  . 

"  =  *  2  ~  kR  ‘  ?  '  H|  (kR  r'°sh  ‘)du. 


(16) 


dR1 


1  e”dkR '  j  y 
?  kTTr~  -  2  / 


H^hkR1  cosh  u)du. 


(17) 
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dR  dR 1 

It  can  be  shown  that  it  and  it  are  not  singular  at  R=0  or  R'=0. 

The  diffracted  electric  field  is  now  in  a  form  which  has  no  singularities 
at  the  source  or  image. 

Since  the  geometrical  optics  terms  have  been  separated  from  the  dif¬ 
fracted  field,  the  impedance  matrix  Z  can,  for  computational  purposes, 
be  broken  into  two  parts;  [Zj  =  [ZQj  +  f AZ  ,  where  ZQ]  represents  the  im¬ 
pedance  matrix  of  a  monopole  on  an  infinite  ground  plane  (E^-0'  contribution) 

[J  i  £ 

AZ  represents  the  E^  field  contribution  which  accounts  for  the  edge. 
Since  PWS  modes  are  used  for  expansion  and  testing  the  ^ZQj  elements  can 

O 

be  computed  in  closed  form  ,  eliminating  the  difficulties  which  would  be 
encountered  if  this  evaluation  were  done  numerically.  This  is  an  important 
and  useful  distinction  from  other  similar  solutions^.  The  remaining  task 
is  to  evaluate  ^Azj. 

The  PWS  modes  are  arranged  on  the  wire  as  shown  in  Figure  2,  numbered 
from  top  to  bottom.  The  test  modes  are  positioned  on  the  wire  axis  at  (x  ,z=0) 
while  the  expansion  modes  are  placed  on  the  wire  surface  at  (x  ,z=a),  where 
a  is  the  wire  radius.  Note  that  the  modes  are  PWS  dipoles  except  for  the 
bottom  mode,  which  is  a  PWS  monopole  with  its  terminal  at  the  base  of  the 
wire.  This  allows  non-zero  current  at  the  attachment  point.  The  contribution 
to  the  mutual  impedance  between  modes  i  and  j  due  to  the  diffracted  field 
can  then  be  written  as 

AZij  =  -  f  W  1  lj(y)Eyjf(y lyo)d*dV  (18) 

yoi  yj 

Thus,  [aZ]  can  be  calculated  with  two  numerical  integrations  for  Equation 
(18)  and  one  for  Equations  (16)  and  (17).  These  integrations  are  well  behaved, 
typically  requiring  only  10-20  Simpson  intervals  each.  This  ease  of  computation 


i 


r  1 

of  M  |  is  a  direct  result  of  the  fact  that  the  incident  and  reflected  field 

C  J 

components  were  separated  from  the  total  field,  allowing  the  nearly  singular 

portion,  fz  j,  to  be  evaluated  in  closed  form. 

The  solution  described  above  gave  a  quickly  converging  result  for  the 

wire  currents,  as  shown  in  Table  I  below.  This  table  qives  AZ  =  Z-  -Z  , 

3  1  n  o 

where  Zjn  is  the  input  impedance  of  a  A/4,  a= .001 A  monopole  near  the  edge 
of  a  half-plane,  and  Z0  is  the  input  impedance  of  the  same  antenna  mounted 
on  an  infinite  ground  plane.  xQ  is  the  distance  of  the  antenna  from  the 
edge.  Results  are  shown  for  1,  2,  and  3  modes  on  the  monopole.  ^  delta 
gap  generator  at  the  base  of  the  monopole  was  used. 


Table  I 


X0/A 

AZ,  1  mode  1 

AZ,  2  modes 

AZ,  3  modes 

.05 

i 

12.l-j7.lQ 

11.8-j8.7Q 

11.7-j9.0Q 

.10 

4.2-J8.2Q 

3.5- j9.2Q 

3.3-j9.4Q 

.20 

-2.8-j3.8Q 

-3.6-j3.8Q 

-3.7-j3.8Q 

.40 

-.32+jl.9Q 

-.09+j2.2Q 

- .04+ j2 . 2Q 

.60  | 

.87-j.730 

.85- j . 91Q 

.85- j.  95.0 

The  convergence  of  ZQ  is  shown  in  Table  II  below: 


Table  II 


Modes 

Zo 

1 

36.5+j21.W 

2 

40.6+ j20.60 

3 

41 .4+j21 .000 

4 

4 1. 7+ j 2 1.30 
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From  the  above  data  it  is  seen  that  the  solution  has  practically  con¬ 


verged  with  two  wire  modes.  Computation  time  for  this  solution  is  about 
30  seconds*  when  2  expansion  modes  are  used.  This  solution  is  valid  for 
any  distance,  x  ,  from  the  half-plane  edge  except  when  the  antenna  is  within 
a  few  wire  radii  of  the  edge.  In  this  case  numerical  problems  are  encountered 
since  the  fields  contain  the  l//x  edge  singularity.  Plots  of  Z-p  vs.  xQ 
are  shown  in  Figures  3  and  4,  where  comparison  is  made  with  results  from 
the  MM/GTD  method  ,  and  with  results  from  the  MM/surface  patch  solution  pre¬ 
sented  in  the  next  section. 
fr- variation  of  wire  current 

The  solution  in  the  previous  section  used  an  axial  current  filament 
for  testing  and  a  circumferentially  uniform  surface  current  located  on  the 
surface  of  the  wire  for  expansion.  Thus,  the  solution  did  not  account  for 
or  give  any  information  about  the  circumferential  variation  of  f*0  current 
around  the  wire  surface.  This  aspect  of  the  solution  was  of  interest  for 
two  reasons. 

First,  it  was  important  to  determine  if  a  full  account  of  the  circum¬ 
ferential  variation  of  the  wire  current,  due  to  the  edge  proximity,  would 
alter  the  input  impedance  as  determined  by  the  previous  solution,  which  did 
not  account  for  any  variation.  Secondly,  it  was  of  interest  to  determine 
the  circumferential  variation  of  the  current  at  the  base  of  the  wire  in  order 
to  build  a  suitable  expansion  mode  for  wires  attached  near  an  edge,  for  use 
in  the  surface  patch  MM  solution. 


♦ATI  times  reported  here  are  for  Datacraft  6024  -  about  8  times  slower 
than  IBM  360/165. 
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The  c l rcumf erent i al  variation  of  the  axial  wire  current  was  taken  into 
account  by  placing  expansion  filaments  evenly  around  the  surface  of  the  wire, 

Q 

i.e.,  by  constructing  a  wire-grid  model  of  the  wire  .  Each  filament  is  an 
independent  unknown,  and  can  also  be  broken  up  into  separate  modes  along 
its  length.  Test  modes  are  placed  alongside  the  expansion  filaments  spaced 
.OOOOOlA  away.  This  procedure  approximates  true  surface  testing  and  expansion. 

figure  5  shows  a  a/4  monopole,  radius  .0Q1X,  placed  .05X  from  the  edge 
of  a  half-plane.  Also  shown  are  the  positions  of  eight  expansion  filaments 
used  in  this  example.  In  this  example  3  modes  were  employed  in  the  longi¬ 
tudinal  direction  for  a  total  of  3x3  =  24  unknowns. 

The  current  at  the  base  of  the  monopole  is  plotted  as  magnitude  and 

phase  in  Figure  5.  As  expected,  the  data  show  symmetry  between  filaments 

?,  3,  4  and  8,  7,  6,  respectively.  The  data  also  show,  as  might  be  expected, 

more  current  flowing  away  from  the  edge  than  toward  it.  If  the  filament 

currents  in  the  above  table  are  added  to  find  the  total  current,  the  "average" 

input  impedance  can  be  calculated  and  is  found  to  be  identical  to  the  value 

given  by  the  simpler  one-filament  solution  of  section  3.  This  interesting 

effect  occurs  for  other  segmentation  schemes  and  distances  from  the  edge, 

and  is  similar  to  an  effect  reported  in  connection  with  closely  spaced  wire 
9 

antennas  .  The  circumf erenti al  variation  of  the  wire  current  becomes  less 
pronounced  as  one  moves  up  the  wire. 

Expressions  for  surface  currents 

In  this  section  simple  closed  form  expressions  will  be  derived  for 
the  surface  current  density  on  a  perfectly-conducting  half-plane  induced 
by  an  inf  ini tes imal  electric  dipole  on  the  surface  of  the  halfplane  and  near 
its  edge.  These  results  will  be  used  for  synthesis  of  an  expansion  mode 
for  a  wire  attachment  near  an  edge  in  the  surface  patch  MM  solution. 
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lJR|.  =  "  =  'MRl.  „  ’  so 


[„  |  -  I„|  -  2  f  H',2,  (kR  cosh  u)du 

KU=0  !  1>=27r  o 


=  2^  ( x ,  z  |  x  i  z  ) 


=0 

o 


where 


K 

4*  ( x ,  z  |  x  , z  )  =  f  h(2^  (kR  cosh  u)du. 


Then  (23),  (24)  become 


s  x  *  5tt7 


2k  ^  «p(x,z|x  z  ) 
0 


sz=WJ  \-fbzz  H(12)(kR1)  -  2k  gf-  *(x,z|x0,z0)| 


fPPo 


We  now  develop  a  small  argument  approximation  for  4>: 
For  small  x. 


Hj2)(x)  =  x/2  +  j  , 


so  (26)  becomes 


J  k^-cosh  jj 

*  J  7  -nit  ft  rnc 


nkR  cosh  u 


1  2 i  1  2  1 

=  k/pp^  cos  7  (<M>0)  ♦  ^  tan’  — 3—  cos  ?  (*-*„). 


For  the  present  case  lt>=lt'o=0,  so  P-*  and  P0  =  xfj.  The  real  part  of  the  approxi¬ 
mation  to  the  function  >p  is  ignored  in  the  remainder  of  the  analysis  since 
it  is  negligible  in  comparison  to  the  imaginary  part. 

The  last  step  is  to  carry  out  the  differentiations  as  indicated  in 
(27), (2 8).  When  this  is  done  there  will  be  present  terms  which  vary  as  l/« 
due  to  the  charges  which  are  located  at  the  ends  of  the  dipole.  Since  these 
terms  vanish  in  a  physical  solution  which  has  continuous  current,  they  are 
omitted  in  the  final  expressions.  So,  for  a  vertical  electric  dipole  of 
moment  I?,  near  the  edge  of  a  halfplane,  the  surface  currents  are 


Observe  that  J  -*0  with  /x  dependence  and  J  +CD  with  l//x  dependence, 

Jo  j  C. 

as  xvQ  as  required  by  edge  conditions.  Also,  J$x  and  J  contain  the  correct 
1/R  singularity  in  the  vicinity  of  the  source.  It  is  also  of  interest  to 
note  that  Equations  (30),  (31)  are  not  frequency  dependent,  i.e.,  they  are 
quasi-static  results.  However,  retaining  the  first  term  in  the  real  part 
or  higher  terms  in  the  imaginary  part  of  (29)  will  introduce  a  frequency 
dependence,  if  desired.  Three-dimensional  plots  of  Equations  ( 30 ) -(31)  are 
shown  in  Figures  6  arid  7. 

IN.  SURFACE  PATCH  MM  ATTACHMENT  MOOES  EOR  WIRE/PLATE 
ATTACHMENTS  NEAR  OR  AT  EDGES 

This  section  deals  with  methods  of  modifying  the  author's  surface  patch 
modeling  program^  to  handle  antennas  mounted  near  or  at  the  edge  of  plates. 

Originally,  the  program  used  an  attachment  mode  which  had  the  following  features 

53 


a)  Enforced  continuity  of  current  from  wire  to  plate. 

b)  Enforced  the  proper  p/p  dependence  for  surface  current  in 
the  vicinity  of  the  attachment  point. 

c)  Had  no  ^-variation  in  surface  current.  (In  this  section 
is  used  to  denote  azimuthal  direction  around  the  attached 
wire.) 

d)  Worked  well  for  wire  attachments  .1\  or  further  from  an 
edge,  and  for  .U<b<.25A. 

This  attachment  mode,  shown  in  Figure  8a,  failed  for  antennas  mounted 
closer  than  .lx  from  the  edge  as  shown  in  Figures  3  and  4  (as  mode  A), 
probably  for  the  following  reasons: 

a)  Oisk  mode  radius  was  limited  to  the  distance  of  the  monopole 
from  the  edge,  thus  reducing  the  area  in  which  the  p/p  singu- 
1 ari ty  was  enforced. 

b)  Disk  mode  did  not  contain  any  (I)- vari  at  ion  in  current  density, 
which  probably  becomes  increasingly  important  as  the  antenna 
gets  closer  to  the  edge. 

Thus,  a  better  attachment  mode  was  needed  to  handle  antennas  mounted  close 
to  or  at  the  edge  of  plates. 

Attac hme rit  Mode  Development 

Four  expansion  modes,  denoted  as  A,B,C,D,  for  wire'plate  junctions 
will  be  described.  All  four  modes  consist  of  a  disk  monopole  which  rests 
on  the  plate,  jnd  i  wire  monopole  which  overlaps  with  the  connecting  wire. 
Sinct?  this  wire  monopole  is  common  to  all  four  attachment  modes,  it  is  not 
described  in  detail  here. 


54 


Mode  A)  -  This  is  the  original  disk  monopole  described  in\  having  a  surface 
current  density 

jp(p,+)  ■  a  ,  (32) 

where  A  is  a  normalization  constant  chosen  to  satisfy  (33): 

?TT 

/  Jp(a,<t>)ad<)>  =  1  amp  .  (33) 

<J)=0 

In  (32),  a  is  the  wire  radius,  and  b  is  the  outer  radius  of  the  disk.  For 
this  mode  the  disk  is  circular,  so  b  is  a  constant  w.r.t.  <t>.  Figure  8a  shows 
how  this  mode  would  be  applied  to  a  wire  attachment  near  an  edge.  Mode  A) 
fails  when  used  near  an  edge  for  the  above  mentioned  reasons. 

Mode  B)  -  The  first  approach  toward  improvement  of  disk  mode  A)  was  to  main¬ 
tain  the  disk  radius  at  about  .2A  where  it  overlapped  ^the  plate,  but  cut 
off  the  disk  where  it  overhung  the  plate;  see  Figure  8b.  The  surface  current 
density  for  this  disk  mode  can  be  written  as 

•yp.f)  *  8  SjC-^W-Cl  ,  (34) 

where  B  is  a  constant  chosen  to  satisfy  (33)  and  b ( 4>)  is  the  outer  radius 
of  the  disk,  now  expressed  as  a  function  of  <1>. 

This  type  of  mode  overcomes  some  of  the  drawbacks  of  the  previous  mode 
in  that  it  has  a  4>- var i ation  of  current  and  yields  a  larger  area  where  the 
l/o  singularity  is  enforced.  The  results  obtained  when  using  this  mode  are 
not  s ign i f icantly  better  than  those  obtained  with  the  original  mode  A). 

This  is  probably  due  to  the  fact  that  the  <t>-dependence  is  not  close  enough 
to  the  exact  dependence.  For  instance,  when  p  is  close  to  a,  the  current 
density  of  mode  B)  is  large  but  does  not  vary  appreciably  with  <J>,  while 


F  iqure  5  iri<J  tqu.it  ion  (30)  show  that  significantly  less  current  is  flowing 
toward  the  edge  than  a way  from  it.  Another  point  is  that  this  mode  does 
not  have  a  uniform  limit  as  the  attachment  point  approaches  the  edge.  One 
would  expect  that  the  current  flowing  toward  the  edge  would  smoothly  decrease 
to  zero  as  the  attachment  moves  to  the  edge.  Mode  3)  does  not  have  this 
behavior.  The  disk  modes  presented  below  do  have  this  desired  asymmetry. 

Mode  C)  -  The  next  step  was  to  construct  a  disk  mode  which  had  less  current 
flowing  toward  the  edge  than  away  from  it.  One  distribution  which  does  this 
is 

Jp(p>t>)  r  C  sin  kb(4>),  (35) 

where  C  is  a  constant  chosen  to  satisfy  (33).  This  current  distribution 
roughly  corresponds  to  that  obtained  from  the  exact  solution  for  the  circum¬ 
ferential  current  variation,  as  shown  in  Figure  5.  This  type  of  mode  has 
the  desired  property  of  uniformly  approaching  a  semicircular  disk  as  the 
attachment  point  moves  to  the  edge  (Figure  8).  Also,  this  mode  reduces  to 
the  standard  attachment  mode,  mode  A),  when  the  distance  between  the  attach¬ 
ment  point  and  edge  becomes  greater  than  the  disk  outer  radius  (v\2A),  since 
b(<j>)  =  constant  for  this  case. 

Th'is  type  of  attachment  mode  gave  the  best  results  for  antennas  at 
or  near  an  edge.  Results  from  the  surface  patch  MM  solution  using  node  C) 
with  a  A/4  monopole  near  the  long  edge  of  a  .4Ax.5A  plate  are  shown  in  Figures 
3  and  1.  The  shift  between  the  surface  patch  results  and  the  exact  results 
for  a  semi- inf  ini te  half-plane  in  these  figures  is  due  to  the  fact  that  a 
finite-sized  plate  had  to  be  used  with  the  former  solution.  As  can  be  seen, 
mode  C)  yields  a  significant  improvement  over  mode  A).  Figure  10  shows 
results  using  mode  C)  when  the  monopole  is  mounted  on  the  edge  of  a  plate 
vs.  angle.  56 


Mode  D)  -  The  next  type  of  disk  mode  which  was  developed  was  one  constructed 
from  an  arbitrary  number  of  wedge-shaped  sectors  of  current,  with  current 
distribution  like  that  of  disk  mode  B)  over  each  sector,  but  with  each  sector 
being  an  independent  unknown.  See  Figure  8d,  which  shows  a  4  sector  disk. 

This  type  of  mode  has  the  advantage  of  allowing  the  MM  solution  to  arrive 
at  currents  (both  magnitude  and  phase)  which  best  satisfy  the  integral  equation 
in  a  least-mean-square  sense.  The  resolution  of  the  variation  in  disk  current 
can  be  made  finer  by  increasing  the  number  of  wedge-shaped  sectors  that  the 
disk  is  broken  into.  When  eight  such  sectors  were  used  the  results  for  input 
impedance  agreed  fairly  well  with  those  obtained  using  mode  C).  Note  that 
this  type  of  attachment  mode  would  work  well  for  wires  attached  to  piecewise- 
curved  surfaces. 

Fdge  Mode  -  Recently,  there  has  been  an  interest  in  directly  enforcing  the 
edge  condition  in  MM  solutions^.  Thus,  it  was  decided  to  construct  a  mode, 
based  on  the  analysis  of  Section  II,  that  would  enforce  the  edge  singularity 
of  the  surface  current  component  parallel  to  the  edge  of  the  half-plane, 
and  study  its  effect  on  the  solution.  This  "edge  mode"  was  constructed  from 
PWS  current  filaments  to  embody  the  essential  features  of  Equation  (31). 

The  surface  current  density  of  this  mode  is 


-sin  k(z-Ht) 
sin  kU-zQr 


sin  kz 


o 


sin  k(z-j) 
sin  kfz0-i) 


for  -l<z<-zQ, 

for  -z  <z<z„ 
o  o 

for  zQ<z<i 


(36) 
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where  the  coordinate  system  is  shown  in  Figure  1,  A  is  a  suitable  normalizing 
constant,  and  2 l  is  the  overall  length  of  the  mode.  The  width  of  the  mode 
(extent  in  x-direction)  is  arbitrary,  but  since  most  of  the  current  is  located 
near  the  singularity  at  x^O,  this  width  is  not  too  critical  if  chosen  to 
be  .05^  or  greater.  The  result  of  using  this  edge  mode  with  different  lengths 
and  widths  was  interesting.  First,  the  edge  mode  was  excited  in  the  surface 
patch  MM  solution,  as  indicated  by  a  non-zero  mode  current.  However,  the 
input  impedance  did  not  significantly  change  from  the  result  obtained  without 
using  an  edge  mode.  Thus,  all  of  the  other  surface  patch,  wire,  and  attach¬ 
ment  currents  changed  in  the  presence  of  the  edge  mode  in  such  a  way  as  to 
leave  the  input  impedance  unaltered.  At  the  present  time  the  best  explan¬ 
ation  we  have  for  this  effect  is  that  the  singular  behavior  of  the  edge 
current  can  be  approximated  adequately  in  a  least-mean-square  sense  by  the 
surf  ace-patch  dipoles  alone.  Incorporation  of  a  special  edge  mode  explicitly 
satisfying  the  edge  condition  does  not  improve  on  this  situation  in  the  present 
case. 

IV.  SOLUTION  FOR  A  RADIAL  MONOPOLE  AT  THE  EDGE  OF  A 

HALF-PLANE  USING  EXACT  GREEN'S  FUNCTION 

Theory 

The  solution  for  a  monopole  mounted  at  the  edge  of  a  half-plane  is 
similar  to  the  solution  described  in  Section  II,  except  that  here  the  exact 
Green's  function  is  expressed  in  spherical  coordinates.  This  alleviates 
the  problem  of  moving  the  vertical  monopole  up  to  the  edge,  while  allowing 
the  monopolo  to  be  oriented  in  a  position  other  than  vertical.  Figure  9 
shows  the  geometry.  The  monopole  lies  in  the  x-y  plane  at  an  angle  si  from 
the  half-plane  surface.  Without  loss  of  generality,  a  is  restricted  to  0«j<ti. 
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From  another  paper  by  Senior  ,  the  radial  component  of  the  electric 
field  due  to  a  radially  directed  electric  dipole,  with  moment  IH,  near  the 
edge  of  a  half-plane  is 


E  .  ,a0  l-jfj  t  2\  ro  , 
Err  '  JiT  |  r  '  ^7  *  )  —  U»'IR' 


where  IR  and  IR,  are  defined  in  (3),  (4).  Equation  (37)  represents  the  fields 
in  terms  of  a  Debye  potential,  whereas  Equation  (2)  represented  the  fields 
in  terms  of  a  Hertz  potential. 

Now,  it  can  be  shown  that 


Lrr  =  7^T< 


;  ,  r°  e-Jf 

.2  *  )  r  TT 


,  and 


Eimage  _  l/Y_  +  k2^  r 

Err  -  4„j  \{dr2  +  kJ  r 


ro  e-J'kR' 


where  the  minus  sign  in  (39)  is  due  to  the  fact  that  the  image  dipole  is 

/V 

in  the  -r  direction.  Then,  defining 


E  =  Finc  +  rima9e  +  Fdi^ 
rr  crr  trr  tr»>  * 


we  have 


where  rrclR  and  irdR  are  defined  in  (16)  and  (17).  Observe  that  when  a<m/2 
there  is  no  reflected  Err  field  for  4»tf/2,  so  (41)  does  not  strictly  repre¬ 
sent  the  diffracted  field  in  the  usual  sense. 
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The  source  and  image  fields  are  then  separated  from  the  total  field. 


and  the 


w 


section  of  the  total  impedance  matrix 


is  computed  in  closed 


form  as  before.  The 


AZ 


elements  are  found  by 


AZij  =  '  i"  /  Ii^r0)  E?rf  (rlr0)Ij(r')dr  dr0.  (42) 

roi  rj 

where  rQj  ranges  over  the  test  mode  and  r-  ranges  over  the  expansion  mode. 

The  differentiations  in  (41)  can  be  eliminated  by  integration  by  parts,  which 

J 

also  removes  one  of  the  integrations  in  (42).  Note  that,  because  of  the 
non-zero  wire  radius,  the  electric  field  along  the  surface  of  the  wire  con¬ 
tains  an  Eq  component  in  addition  to  the  Er  component.  This  E^  component, 
which  makes  only  a  small  contribution,  has  been  ignored  here. 

This  solution  was  programmed  and  run  for  a  x/4,  a=.001\  monopole,  with  I 

a  delta  gap  generator  at  the  base  of  the  monopole,  vs  angle  a.  Figure  10 

shows  Z .  using  two  modes  on  the  wire  compared  with  the  surface  patch  VM  i 

solution.  This  result  at  a=iT/2  also  agrees  with  the  analysis  for  the  vertical 

monopole  in  Section  II  as  xQ  approaches  zero.  Note  that  Z-n  goes  to  zero 

as  a  approaches  zero  as  expected,  since  the  antenna  then  becomes  shorted 

by  the  ground  plane.  As  a  approaches  zero,  the  reactance  curve  takes  a  sharp 

turn  downward  before  going  to  zero.  This  effect  is  probably  explained  by 

the  fact  that  the  capacitance  between  the  monopole  and  half-plane  increases 

sharply  as  the  monopole  gets  closer  to  the  half-plane,  before  going  to  zero 

as  it  touches  the  half-plane. 

Another  independent  check  is  available  from  Tai5  who  computed  the  input  -r  J 

’•esi  stance  of  a  monopole  at  by  integrating  the  far  field  radiation  pattern 

(the  input  reactance  cannot  be  calculated  using  Tai's  method).  Also  shown 

60 


I 

S 


•  nmn.  m  jt.  r 


in  Figure  10  is  a  surface  patch  result  obtained  by  Glisson^  for  a  \/4  monopol 
Mounted  on  the  edge  of  a  finite  plate,  which  compares  reasonably  with  the 
other  results. 

The  surface  current  density  can  be  derived  for  a  source  at  the  edge 
of  a  half-plane  in  the  same  manner  as  was  done  in  Section  II.  The  final 
result  for  a  dipole  of  moment  I/gr  is 


l 

r 


.  *0 

sin  y~ 


A/m  . 


(43) 


The  dipole  moment  must  be  reduced  as  /r^  as  r  -*0  to  maintain  finite  fields. 
Note  the  correct  source  and  edge  singularities. 

V.  CONCLUSION 

This  paper  has  presented  analysis  and  results  of  the  problem  of  a 
monopol e  antenna  mounted  near  or  at  the  edge  of  a  half-plane  or  plate.  The 
canonical  solutions  using  exact  Green's  functions  have  been  presented  as 
a  basis  for  comparison  with  other  solutions.  Expressions  for  near-zone 
surface  current  density  were  derived  and  used  to  synthesize  an  attachment 
mode  to  handle  wire  attachments  near  plate  edges  in  a  general  surface  patch 
MM  solution. 
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1.  Geometry  and  coordinate  system  for  a  vertical 
monopole  near  the  edge  )f  a  half-plane. 


Figure  2.  PWS  expansion  mode  arrangement  on  wir 


Input  resistance  of  a  A/4  monopole  vs.  distance,  d,  from  the  edge  of  a  hai*- 
computed  using  MM/exact  Green's  function,  MM/GTD,  surface  patch  with  old  att 
mode  (A),  and  surface  patch  with  new  att.  mode  (C).  ~ 


INPUT  IMPEDANCE  OF  MONOPOLE  (X/4) 
NEAR  EDGE  OF  HALF  -  PLANE 


mndf  (A),  and  surface  patch  with  new  att.  mode  (C). 
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Figure  5.  Magnitude  and  phase  of  filament  current  at  the  base  of 
a  monopole  .05A  from  a  half-plane  edge,  vs.  position 
around  the  wire.  The  layout  of  the  eight  expansion 
filaments  is  also  shown. 
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Figure  7.  Three-dimensional  plot  of  the  magnitude  of  the  surface-current 
density  (Equation  (29))  for  a  vertical  infinitesimal 
electric  dipole  .05A  from  the  edge  of  a  half-plane. 
.00lA<x<.lA,-.  125a<z<.  125A. 


Kijinal  attachment  mode,  mode  A) 


attachment  mode  C)  when  monopole  is 
nounted  at  a  plate  edge. 


final  attachment  mode,  mode  C) 


(d  ) 


segmented  attachment  mode,  mode  D), 
shown  with  4  segments. 


Figure  9.  Geometry  and  coordinate  system  for  a  monopole  mounted 
at  the  edge  of  a  half-plane. 
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